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I. INTRODUCTION 

The discovery of neutrino masses and mixings through solar , atmospheric |3[, and 
reactor J] neutrino oscillation experiments has provided new clues to solving the mysteries 
of the Standard Model. Flavor mixing in the lepton-sector, together with the well-known 
mixing in the quark-sector, may lead to a new understanding of what 'flavor' is, why there 
are three generations of fermions, and where the quark- and lepton-mass hierarchies come 
from. CP violation in the neutrino-sector could potentially be large enough to account for 
the matter-antimatter asymmetry in the universe j^. 

Due to these possibilities, several long-baseline (LBL) neutrino oscillation experiments 
are in progress p, Q, 0, 0, with many more being proposed n[ 12 , 3, 14, 15, 16. Il7| 
for the purpose of better determining the neutrino masses and mixing parameters. Since 
the neutrino beams of LBL experiments necessarily traverse the Earth, the understanding 
of matter effects Q| is crucial in extracting the masses and mixing parameters in vacuum 
from the measured oscillation probabilities in matter. 

The calculation of matter effects requires the diagonalization of the effective Hamiltonian 
in matter, which is an energy- and matter-density-dependent 3x3 matrix. This can be done 
numerically on a computer with ease, but the dependence of the oscillation probabilities on 
the vacuum parameters is not transparent Q|. It is also possible to write down exact 
analytical expressions for the effective parameters in matter in terms of those in vacuum 
[20) 1 . but the expressions are too complicated to be illuminating. Approximate expressions 
for the oscillati on p robabilities have also been worked out, but were limited in their range 
of applicability j^lj], or still too complicated to be of practical use 

In this paper, we derive simple approximate expressions for the effective mass-squared- 
differences and mixing angles in matter, in terms of the corresponding parameters in vacuum. 
They are simple enough to be used by hand, yet accurate enough in a wide energy range. 

This paper is organized as follows. In section II, we review the formalism of neutrino 
oscillations to fix our notation, and list the current experimental bounds on the mixing 
angles and mass-squared differences. In section III, we derive our approximate expressions 
for the effective mixing angles and mass-squared differences in matter, for both neutrinos 
and anti-neutrinos. In section IV, we provide sample calculations to illustrate the accuracy 
of our approximation. In section V, we present simple analytical approximations for the 



oscillation probabilities and demonstrate their utility in understanding which parameters 
are best constrained at what baseline lengths. Section VI concludes. 



II. BASICS 



The MNS Matrix 



Assuming three-generation neutrino mixing, the flavor eigenstates \u a ) (a = e,/i, r) are 
related to the three mass eigenstates (i = 1, 2, 3) via the Maki-Nakagawa-Sakata (MNS) 
matrix: j3| 

(VMNs)ai = (v a \Vi) , \n) = ^2 (^MNs)m Wa) , Wa) = ^ (^MNs)aj Wi) ■ (1) 

a=e,/i,T j=l,2,3 

A popular parametrization is given by 



with 
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(2) 



(3) 



where s y - = sin0 y -, and Cy = cos Without loss of generality, we can adopt the convention 
< Oij < 7r/2, < 5 < 2tt (2^. Of the six parameters in this expression and the three 
masses, which add up to a total of nine parameters, neutrino— ^-neutrino oscillations are only 
sensitive to six: 

• the three mixing angles: #12, #23, #13, 

• two mass-squared differences: Sm^i, 5m| 1; where fimjj = mf — m 2 , and 



the CP-violating phase: 5. 



The Majorana phases, a,\ and 02, only appear in lepton-number violating processes such as 
neutrinoless double beta decay, and cannot be determined via neutrino^neutrino oscilla- 
tions. The absolute scale of the neutrino masses also remain undetermined since neutrino 
oscillation is an interference effect. 



B. Oscillations 



If a neutrino of flavor a is created at x = with energy E, then the state of the neutrino 
at x = is 

3 

Wa,o(x = 0)) = \v a ) = ^(Fmns)^ Wi) ■ (4) 



8=1 



At x = L, the same state is 



\V« )0 (X = L))=J2 (^MNs)aj Wi) ■ 



i=l 



Using 



we find 



Pi = JE 2 -m* = E 



2E 



+ ■■• , 



Wa,o{ X = L )) = e 



EL 



X^exp (-^l) (Vmns) 



Wi) 



(5) 



(6) 



(7) 



Therefore, the amplitude of observing the neutrino of flavor (3 at x = L is given by (dropping 
the irrelevant overall phase) 



{vpWa,o{x = L)) 



3=1 
3 



^(Vmns)/Jj exp f -i^L J (Vmns~ 
7=1 ^ ' 



aj > 



(8) 



and the probability of oscillation from \v a ) to |z^) with neutrino energy E and baseline L is 
given by 

2 



y exp ( -i-^gL ) (V MN 



<W - 4 £ ^(U^UptUaU^) sin 2 ^ + 2 £ ^(U^U^) sin A 



A, 



(9) 

where 

<5m?- Smf; (eV 2 ) „ „ 9 

Ay = -^L = 2.534 E ^ v) ; L (km) , 5m 2 . = m 2 - m) . (10) 

Since 

A 32 = A 31 - A 21 , (11) 

only two of the three Ay's in Eq. JSJ) are independent. Eliminating A 32 from Eq. Q, we 
obtain 

i>« - *„) = 1 - 4 |f/ Q2 | 2 (1 - |£/ Q2 | 2 ) sin 2 ^ - 4 |[/ a3 | 2 (1 - |£/ a3 | 2 ) sin 2 ^ 

+ 2|^ Q2 | 2 |f/ a3 | 2 (4sin 2 4r sm2 ^r + sinA 2i sin A 3 i J , (12) 



for the a — (3 case, and 



2„- 2 A21 . , , rr , 2 i r r |2 • 2 A 3 1 



%^^) = 4|^ 2 | 2 |^ 2 | 2 sin 2 ^ + 4|^ 3 | 2 |[/ / 



/33| 



sin 



2 



+2 RiU^UpaUrtUfc) \4 sin 2 ^ sin 2 ^ + sin A 21 sin A 31 
+4 J( Q ,/3) ^sin 2 sin A 31 - sin 2 sin A 2i ^j , (13) 
for the a 7^ (3 case, where e/( a ,^) is the Jarskog invariant: 

^ = +%(K 1 u f31 u a2 u; 2 ) = +%(u* a2 u f32 u a3 u; 3 ) = +^(u* a3 u /33 u al u; 1 ) 

= -QiT&UfnUriUh) = -Smv^UfnUrtUfc) = -Z(U* a3 U p3 U a2 u; 2 ) 

= ~J(P,a) ■ (14) 

In the parametrization given in Eq. (J5J), we have 

J(ji,e) = —J(e,n) = J(e,r) = ~ J(r,e) = J(t,h) = ~J(/i,t) = A Sm 5 , (15) 

with 

A = si 2 ci 2 si 3 c 2 3 s 23 c 23 . (16) 

The oscillation probabilities for the anti-neutrinos are obtained by replacing U a i with its 
complex conjugate, which only amounts to flipping the sign of 5 in the parametrization of 
Eq. (JHJ). It is clear from Eq. ((TJj) that P{p a v a ) = P(u a — > is a ), which is to be expected 
from the CPT theorem. For flavor changing oscillations, only the Jarskog term in Eq. (|13|) 
changes sign. 



C. Current Experimental Bounds 



When |A 2 i| <C |A 31 | = 0(1), the expressions given above can be expanded in A 21 to 
yield 
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where 

|^ 2 | 2 g (UZzUfaUcaUh) f . 

Kaa ~ i - \u a3 \ 2 ' Ka "~ \u a3 \ 2 \u P3 \ 2 ■ [U) 

Neglecting terms of order 0(A| X ) and higher, we obtain the following simplified expressions 
for a few specific processes: 

P(Ve - Ve) = 1 - 44(1 - 4) Sin 2 ^ " 

= 1 - sin 2 (20 rct ) sin 5 
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31 — ^12^21 



A,i - s? A 



P(u, - v,) = 1 - 4c? 3 4 3 (l - «) ^ 2 ( A31 2 Kw,Aai 

= l-sin 2 (2g atm )sin 2 ( A31 ~ 2 ^ A21 ) , 

v e) = 4(C 13 S 13 S 2 3 - J( / x, e )A 2 l) Sill I 

= 4 (sin 2 6> rct sin 2 6> atm - A sin 5 A 2 i) sin 2 ^— — 



1 j 
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where 
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Kfjte 


s 2 _ 



\ ILL 

We have made the identifications 



sin(2^i 2 ) sin(20 rct ) sin(26> atm ) y/l - tan 2 # rct tan 2 9 atm , 

: 2 2 - s 2 12 ) tan 2 # rct tan 2 # atm - ( -— 2 j cos 5 , 

\ COS t7 rct COS (7 atm / 

~ TTi ^75 )cos5. (21) 

sin 9 Ict sm 9 atm J 



sin 6> atm = s 23 ci3 = sin 9 23 cos 13 , 

sin 9 Ict = S13 = sin 9 13 , (22) 

where 9 at m and 9 TCt are the mixing angles extracted from atmospheric and reactor jj] 
neutrino oscillation experiments, respectively, based on two-flavor oscillation analyses. The 
current experimental bounds from atmospheric neutrinos at the 90% confidence level are 

\5m 2 31 \ = (1.5 ~ 3.4) x 10" 3 eV 2 , 
sin 2 (20 atm ) > 0.92 . (23) 

Only the absolute value of 8m 31 is known since Piy^ — > v^) = 1 — sin 2 (2# atm ) sin 2 (A3i/2) at 
leading order in A21. The 90% confidence limit on 9 rct = #13 from the CHOOZ experiment, 
which measured P{p e — > u e ) = 1 — sin 2 (26 l rct ) sin 2 (A3i/2), depends on the not-yet-well-known 
value of I I : 



sin 2 (2# rci ) < 0.20 for ^m 2 ^ = 2.0 x 10~ 3 eV 
sin 2 (2# rci ) < 0.16 for \5m 2 31 \ = 2.5 x 10~ 3 eV 2 



2 



sin 2 (2fl rci ) < 0.14 for \5m 2 31 \ = 3.0 x 10" 3 eV 2 . (24) 
When |A 31 | > |A 2 i| = 0(1), P{v e — > u e ) simplifies to 

P(v e ^is e ) = 1 - 2|f/ e3 | 2 (1 - |f/ e3 r) -4|f/ el | 2 |f/ e2 | 2 sin s 



|2 • 2 ^21 



\8ml\ (eV 2 ) 


/ 3 m 21 


Upper bound on #13 


2.0 x 10" 3 


0.20 ~ 0.21 


0.23 


2.5 x 10~ 3 


0.18 ~ 0.19 


0.21 


3.0 x 10~ 3 


0.16 ~ 0.17 


0.19 



TABLE I: Comparison of the size of the ratio Jm^/I^m^ | and the 90% confidence limit on #13. 
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(25) 



Since S13 is known to be small from CHOOZ, we can identify 9 12 with the 2-flavor solar 
mixing angle 9^. The current experimental bounds at the 90% confidence level from solar 
neutrinos are p, Q] 



5m 2 21 = 8.2t°A x 10" 5 eV 2 



-0.5 



tan 2 9 sol 



0.4012:8? . 



(26) 



The sign of Sm 2 ^ = rn 2 , — m\ is known to be positive since mf. > m\ is required for the MSW 
effect II8SI to work. 



D. The Sizes of 6>i 3 , 9 12 , and 9 



23 



If we allow both 5m 2 ± and l^m^l to move within their respective 90% confidence limits 
given in Eqs. (J23|) and (J26|) . the ratio of the two is in the range 



Stti 2 

— - f- = 0.023 ~ 0.059 . (27) 

Thus, the approximation A 2 i <C | A. 31 1 that was used above is justified. The square-root of 
this ratio, which we will call e, is in the range 



5m 2 21 



|fo4| - 0.15 - 0.24. (28) 

Note that e is roughly in the same range as that of the 90% upper bound on #13 which can 
be obtained from Eq. (J24j) . In fact, the two numbers are positively correlated as shown in 
TABLE I. 



While we currently have no experimental lower bound on # 13 , many analyses of LBL 
experiments only consider cases where sin 2 (26*i3) > 0.04, which corresponds to #13 > 0.1, 
since smaller values of # 13 would put Piv^ — > v e ) below detectable range. Therefore, we will 
assume that 6^3 is of order e in the following. However, the formulae we derive below can 
be applied as is to cases in which #13 is smaller, since we only use the size of #13 to decide 
when terms containing S13 can be neglected. A smaller 13 will simply make those terms 
even more negligible. 

The 90% confidence limits on the solar mixing angle given in Eq. ()26|) translates into 

(0.18±0.01)vr , (29) 

0.50 ~ 0.57 , 
0.82 ~ 0.87 , 
0.86 ~ 0.94 , 

0.34 ~ 0.50 . (30) 

Being sines and cosines, the upper end of these ranges are always smaller than one. However, 
S12, C12, and sin(2#i2) are still much larger than e so we will treat them, and also #12, as 
numbers of order 1. On the other hand, cos(26 , 12 ) is only slightly larger than e. Its range is 
roughly equal to that of 2e. Therefore, we can consider cos(26'i2)/2 as a number of order e. 
The 90% confidence limits on the atmospheric mixing angle given in Eq. (|2*3*|) translates 

to 

9 atm = 0.64 ~ 0.93 = (0.25 ± 0.05)vr . (31) 

Though we made the identification sin # atm = S23C13 in Eq. (j22j) . the limits on 6*23 are virtually 
identical to those of # atm due to the smallness of #13. We can therefore assume 

s 23 = 0.60 ~ 0.80 , 
c 23 = 0.60 ~ 0.80 , 

sin(2#2 3 ) > 0.96 , 
I cos(26> 23 )| < 0.28 . (32) 

As in the case of 6*12, we can assume s 2 3, c 2 3, sin(26 , 2 3), and 6*23 to be numbers of order 1, 
while cos(2#23) is of order e or smaller. 



e 12 = o.56i£S = 

and 

S12 = 
C12 = 
sin(2fl 12 ) = 
cos(2# 12 ) = 



III. MATTER EFFECTS 



A. Diagonalization of the Effecive Neutrino Hamiltonian 



If the matter density along the baseline is constant, matter effects on neutrino oscillations 
can be taken into account by simply replacing the MNS matrix elements and mass-squared 
differences with their "effective" values in matter: 
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IJ 5 



(33) 



where U is the unitary matrix that diagonalizes the modified Hamiltonian, 
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(34) 



and 



Xi - Xj . 



(35) 



2E~ 7 ~ v 

The factor a is due to the interaction of the \u e ) component of the neutrinos with the 
electrons in matter via ly-exchange: 



2V2G F N e E = 7.63 x I(T 5 (eV) 2 ( ( j^-) ■ 

yg/crcry \GeV / 



(36) 



Note that a is -E-dependent, which means that both U and A.^ are also -E-dependent. It 
is also assumed that E My/ since the ^-exchange interaction is approximated by a 
point-like four-fermion interaction in deriving this expression. 
To see the effect of a, we introduce the matrix 



Q = diag(l,l,e w ) , 
and start with the partially diagonalized Hamiltonian: 
H' 



(37) 
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C12S12C 



13 



s 2 r 2 



C12C13S13 S12C13S13 



Sl2 c 13 s 13 



'13 



aC 12 C 13 



OCi2Si2C 13 OCi2Ci3Si3 

aci 2 si2C 2 3 as 2 2 c 2 3 + 5m 2 2l asi 2 Ci 3 Si 3 
aci 2 ci 3 si 3 as 12 ci 3 si 3 as\ 3 + 5m\ x 



(38) 



The matrix Q serves to rid H' of any reference to the CP violating phase 5. Our strategy 
is to approximately diagonalize H' through the Jacobi method using e = \J ~5rr^J\5rr^\ as 
the expansion parameter. We assume 6 13 = 0(e) as discussed above. Corrections to the 
eigenvalues of H' of order £ 3 |5m| 1 | and higher, and those to the elements of U of order e 3 
and higher will be neglected. For e = 0.15 ~ 0.24, we have e 3 = 0.0034 ~ 0.014. 
Recall that for 2x2 real symmetric matrices, such as 



M 



a j3 
(3 7 



(39) 



diagonalization is trivial. Just define 



R 



where 



cosu; 



tan 2u = 



2/3 
7 — a 



(40) 



and we obtain 



R ] MR 



with 



Ai 



IK 



c 2 , - s 2 , 



Ai 
A 2 

A2 



as: 



(41) 



(42) 



The Jacobi method entails iteratively diagonalizing 2x2 submatrices of a larger matrix in 
the order that requires the largest rotation angle at each step. In the case of H', as we will 
show below, two iterations at most are sufficient to achieve approximate diagonalization to 
the order required, regardless of the size of a. 



In the following, we will consider the five cases a/l^mU = 0(e 3 ), 0(e 2 ), 0(e), 0(1), and 
0(e~ l ) separately. Values of a/l^ra^l outside this range can be treated in a similar manner 
as the 0(e 3 ) and 0(e~ l ) cases. (For conditions on the Earth, however, a/l^mU > 0(e~ 2 ) 
implies E > M\y, invalidating our approximation.) 



Al. a/ IJml! | = 0(e 3 ) case 



We begin by considering the case a/l^mjjj = 0(e 3 ). Treating both c\ 2 and s\ 2 as numbers 
of order 1, and s 13 as a number of order e, the relative sizes of the elements of H' are given 
by 



H' 



ac 12 c 13 
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ac 12 s 12 c l3 as 12 C 13 ^ m 21 as 12 c 13 s 13 
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31 1 



0{e 3 ) 0{e 3 ) 0{e 4 ) 
0(e 3 ) 0(e 2 ) 0(e 4 ) 
0(e 4 ) 0(e 4 ) 0(1) 



(43) 



OCi2Ci 3 Si3 CtSi2Ci3Si3 <2S 13 

Of the three 2x2 submatrices of H ; , the one which is farthest from diagonal is the (1, 2) 
sub-matrix since it requires an 0(e) rotation to diagonalize, while the other two submatrices 
only require 0(e 4 ) rotations. To diagonalize the (1,2) submatrix of if', we define 
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1 



(44) 



where 

c v = cos <p , 
Using V, we find 



sin ip , tan 2cp = 



ac 2 n sin 29i 2 



8m 2 i — ac 2 3 cos 29 



12 



< f < 
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(45) 



H" = V^H'V 



K 






An 
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asi 2 ci 3 si 3 
ac' 12 c 13 s 13 asi 2 ci 3 si 3 as\ 3 + Sm^ 
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where 
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sin# 



12 ' 



e' 12 = e 12 + ip, 



(47) 



and 
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(as 2 2 c 2 3 + Sm 2 21 )cl - (ac 2 2 c 2 3 )sl 
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+ Sm 2 ,-^ ± a/ (acf 3 — b~m 21 ) 2 + 4acf 3 sf 2 5m 



21 



tan 26C 



tan 2^12 + tan 2ip 



Sm?,! sin 2^12 



12 1 — tan 2# 12 tan 2cp 5m 21 cos 2# 12 — ac\ 3 ' 
which can be used to calculate 9' 12 without calculating (p first. 



(48) 

(49) 
(50) 



Since a/\8m\i\ = 0(e 3 ), which means that a/Sm^i = 0(e), we can expand A' 1; \' 2 , and <p 
in powers of a/5m 21 and s± 3 = 0(e). We find 

A; = ac\ 2 + 0(e i \5m 2 3l \) = 0(e 3 \5m 2 1 \) , 



A' 2 = 5m 2 2l + as\ 2 + 0(e i \5m 2 3l \) = 6m 2 21 + 0(e 3 \5m 2 31 \) = 0(e 2 \5m 2 31 \) 



(51) 



and, 
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13 



5m 2 21 



sin(2# 



12) 
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13 
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sin(2fl 12 ) cos(2fl 12 ) + ■ ■ 



2 



sin(2# 12 ) + 0(e 3 ) = 0(e) , 



12 — Q\2 + <fi — $12 + 



2 8m 21 



sm(29 12 ) + 0(e 3 ) = 0(1) . 



(52) 



Note that the second term in the expansion of <p can be considered to be of order e 3 since 
cos(26 l i2)/2 is of order e as discussed in the previous section. Therefore, the sizes of the 
elements of H" are 



H" 



ac 12 ci 3 s 13 
j 



2 I 
31 1 



0(e 3 ) 0(e 4 ) 
0(e 2 ) 0(e 4 ) 
0(e 4 ) 0(e 4 ) 0(1) 



(53) 



a; o 

A 2 as' 12 ci 3 si3 = \Sm 

ac' 12 c 13 s 13 as' 12 c 13 s 13 as 2 3 + 5m 2 3A 

The rotation angles required to diagonalize the (1,3) or (2,3) sub-matrices are of order e A , 
which we will neglect. So in this case, the Hamiltonian is approximately diagonalized by 
one (1,2) rotation of angle 0(e), and the eigenvalues are 

Ax « a; 



ac 12 , 



A 2 ~ A' 2 ~ 5m 21 + as 



2 

12 ) 



A, 



\ 3 « 5m 2 31 . (54) 

We have kept terms up to O^IJmjj-J) here to show how the eigenvalues are shifted away 
from their vacuum values. 



A2. a/\5m%y\ = 0{e 2 ) case 



Next, we consider the case a/lSm 2 ^ 
this case are 



0(e 2 ). The relative sizes of the elements of H' in 
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0(e 2 ) 0(e 2 ) 0(e 3 ) 
0(e 2 ) 0(e 2 ) 0(e 3 ) 
0(e 3 ) 0(e 3 ) 0(1) 



(55) 



and again we find that we must diagonalize the (1,2) submatrix first. We define V and ip as 
in Eqs. ()44|) and (J45|) . and the matrix will be partially diagonalized to H" in Eq. ()46|) . with 
X[, A' 2 , and 9' 12 defined as in Eqs. (147)1 and (|48|1 . Since a and Sm 2l are of the same order in 
this case, both tan(2<£>) of Eq. ()43j) and tan(26^ 2 ) of Eq. (joTIj) can be expected to be large. 
Therefore, both s' 12 and c' 12 can be considered to be numbers of order 1 in this case also, 
and the relative sizes of the elements of H" are 



H" 



a; 

A' 2 as' 12 ci 3 si3 

ac' 12 c 13 s 13 as' 12 c 13 s n asj 3 + Sm^ 



ac l2 c 13 s 13 
j 



\8m 



2 I 

31 1 



0{e 2 ) 0{e 3 ) 
0{e 2 ) 0{e 3 ) 
0(e 3 ) 0(e 3 ) 0(1) 



(56) 



Further diagonalization requires rotations by angles of order e 3 , which we will neglect. 

In this case, we cannot expand ip, X'±, and A' 2 in powers of a/Sm 21 or its inverse. However, 
we can still expand in S13 = 0(e) and find 

(a + 5m 2l ) ± yjja — 5m 21 ) 2 + 4a £m 21 s 2 



X' 



'12 



+ 0(6*1^1) 



and 



<P 



tan 



a sin 29 



12 



8m 



21 



a cos 29 



12 



a <5m 21 s 12 ci 2 s 2 3 



[a - 5m 2 21 ) 2 + 4a Sm^s 2 + ° ^ ' 



'12 



The coefficient of s\ z in the second term is bounded from above by 



a5m 21 si 2 ci2 



< 



C12 



0.38 ~ 0.42 » 2e 



(57) 



(58) 



(59) 



(a — <5m 21 ) 2 + 4a 5m 21 sf 2 4s i2 
Therefore, though this factor is formally of 0(1), we can consider it to be a number of order 



e and approximate 



if ~ - tan 1 
^ 2 



a sin 2#i 2 



<5m 



21 



a cos 2$ 



12 



^^tan" 1 f/ 2 mllSin f 12 
12 2 V^2i c os2# 12 -a 



(60) 



The eigenvalues in this case are 



a 2 « a; 



+ (Jmli) — \/ (a — Srn^) 2 + 4a Srrv^s 



2 

12 



a + Sm^x) + a/(o — ^m^) 2 + 4a Jm^s^ 



A 3 « 5m 2 3l . (61) 



y!5. o/l^mgil = 0(e) case 

The relative sizes of the elements of H' in this case are 



H' 



ac 12 c 13 



OCi2Si2C 13 CtCi2Ci3Si3 



aci 2 Si 2 c 2 3 as\ 2 c\ 3 + 8m 2 x asi 2 Ci 3 Si 3 



aci 2 ci 3 si 3 asi 2 ci 3 si 3 as? 3 + <5m 



'31 



\5m 2 31 \ 



0(e) 0(e) 0(e 2 ) 
0(e) 0(e) 0(e 2 ) 
0(e 2 ) 0(e 2 ) 0(1) 



• (62) 



As in the previous two cases, the (1, 2) submatrix is diagonalized by V defined as in Eq. (I44J) 
with tp defined as in (145)1 . The partially diagonalized form is H" in Eq. (146)1 . with A' 1; A' 2 , 
and 9' 12 defined as in Eqs. (gZJ) and (gBjl - 

Since Sm^/a = 0(e) in this case, we can expand A' 1; A' 2 , and in powers of Sm 2l /a and 
find 



Ai 

X'o 



2 - 2 + 0(e 3 5m 2 31 ) 



5m 21 c 12 



5m 01 s 



2 

21°12 



0(e 3 Sm 2 31 



0(e 2 \6m 2 31 \) , 
= 0(e\5m 2 31 \) 



(63) 



and 



<P 



#12 
#12 



9' 



'12 - #12 + V? 

which shows that 



12 



'12 



ac 



1 / 5m 2 ,! 

2 ' 

^21 

2a 



13 



sin(2# 



12 , 



sin(2# 12 ) + 0(e 5 



7l 



5m. 



21 



cos — — 



sin 



2 " 


2a 


71 


5m 2 2l 


2 ~ 


2a 


7T 




2 ~ 


2a 



sin(2# 12 ) + 0(e 3 



sin(20 12 ) 
sin(2# 12 ) 



ac 



13 



sin 



sin(2# 12 ) cos(2fl 12 ) + 



2a 



sin(2fl 



12 



I sm(26'i2, 



0(e), 
0(1). 



(64) 



(65) 



Note that 



ac 



12 



5m 2 2l 



sin(2# 



12 



1 5m 



2 I 
31 1 



5m 



sin(2fl 12 ) = |5tt4|0(£ 2 



2a '"'J \2\5m 2 n \ 

for all a 3> 5m 2 i- We will use this relation repeatedly in the following. 
Thus, we find the sizes of the elements of H" in this case to be 



(66) 



H" 



K 
o 





a; 



ac 12 C 13 s 13 

as' 12 ci 3 si3 
ac' 12 c 13 s 13 asi 2 ci 3 si 3 as 2 13 + 5m\ l 



\5m 



2 I 
31 1 



0{e 2 ) 0(e 3 ) 
0(e) 0(e 2 ) 
0(e 3 ) 0(e 2 ) 0(1) 



(67) 



This time, the (2,3) submatrix requires an 0(e 2 ) rotation to be diagonalized. We define 



W 



where 



= cos i 



smi 



1 

S0 

—Sf c<f, 

, tan 20 \ 



(68) 



as' l2 sin 26 



13 



5m 2 31 



as 



13 



A' 



(69) 



The angle is in the first quadrant when 5m 2 1 > 0, and in the fourth quadrant when 
5m\i < 0. Then, 



H'" = W ] H"W 



a; 

-ad 12 c 13 s 13 s^ 

ac/ 12 c 13 s 13 c </> 



-ac'^cns^s^ ac' 12 ci 3 si 3 C0 



o 



A 3 



(70) 



where 



AS 



A 3 



31; 



r 2 - s 2 



(as 2 3 + 5m 2 1 )c 2 - \' 2 s 2 



"31 , 



c 2 - s 2 



(71) 



If we define 



A ± 



[A ' 2+(as 2 3 + 5m 2 i)]± /[A , 



as 



13 



5m 2 1 )] 2 + 4a 2 S ' 12 2 c 2 3S 2 3 



(72) 



then 



A 2 = A 7 ! , A 3 ' = A'j , if (5m 2 ! > , 



K — A+ ; A3 — A'_ 



if Sm^ < . 



(73) 



Expanding A 2 , A3, and in powers of a/Sm 31 = 0(e), we find 



X'2 = \' 2 + 0(e 4 \5ml 1 \) = a + 5m~ 21 s 12 



(5ml, + as 2 13 ) + 0(£ 4 |5m 2 31 |) 



+ 0(e 3 |^ 1 |) = 0(^1^1), 
0(e 3 |^|) = 0(\8m 2 31 \) 



5m 2 31 



(74) 



and 



sfn(2# 13 ) + 



5m 2 3l 



'13 



0( £ 3 ) = 0(e 5 



which means that = 0(e 1 



0(1), and 



H" 



Ai -ac'^cississ^ aci 2 ci 3 si 3 c^ 



-ac 12 ci 3 si 3 s^ 

ac 12 C 13 s 13 c <?i 







I bra 



(75) 



0(£ 2 ) 0(£ 5 ) 0(£ 3 ) 

0(e 5 ) 0(e) 
0(e 3 ) 0(1) 



(76) 

Further diagonalization require rotations of order 0(e 3 ) and higher, which we neglect. So in 
this case, the Hamiltonian is approximately diagonalized by a (1,2) rotation of angle 0(1) 
followed by a (2, 3) rotation of angle 0(e 2 ), and the eigenvalues are 



Ax « a; « 

A 2 ~ A 2 ~ 
A 3 « A'' « 



(5?77, 21 C 12 , 

a + 5m 21 s 2 2 , 
Sm 2 31 . 



(77) 



>L£. a/\5m 2 31 \ = 0(1) case 



The relative sizes of the elements of H' in this case are 



H' 



ac 12 c 13 



aci 2 s 12 c\ 3 aci 2 c 13 s 13 

2 „ n 2 „2 , 



\5m 2 31 \ 



0(1) 0(1) 0(e) 
0(1) 0(1) 0(e) 
0(e) 0(e) 0(1) 



(78) 



aci 2 Si 2 c 13 as 12 c 13 + 5m 21 asi 2 c 13 s 13 
aci 2 ci 3 si 3 asi 2 ci 3 si 3 as 2 3 + Sm^ 

As in the three previous cases, the (1,2) submatrix is diagonalized by V defined as in Eq. ()""""")) 
with (f defined as in (|45|1. The partially diagonalized form is H" in Eq. (|46|). with A' 1; A 2 , 
and 6^ 2 defined as in Eqs. 14*7)1 and (f4*5)l . The expansions of A' x and A' 2 in powers of Sm^/a 
yield 



A', = ^c^ + O^I^I) = 0(e 2 \5m 2 3l \) 



A' 2 = ac? 3 + ^ lS ? 2 + 0( £ 4 |^3il) = 0(\Sm 2 31 \). 



(79) 



Note that we cannot replace C13 in this expression with 1 without introducing an error of 
order e^JmjjJ which is the same order as the second term. The expansions of ip and 8' 12 are 



6' 



12 



7T 

2 

#12 ■ 



" #12 



2a 



sin(2# 12 ) +0(s 4 ) 



|-^sin(2# 12 )+0(. 4 ) 



(80) 



From Eqs. and (|p^j). we can tell that s' 12 = 0(1), and 



12 



0(e 2 \8ml 1 \). Therefore, 



H" 



K 






A 2 



ac 12 C 13 s 13 

as' 12 ci 3 si3 



/ / 2 

ac 12 C 13 s 13 a "S 12 Ci3Si3 as 13 



\8m 



0(e 2 ) 0(e 3 ) 
0(1) 0(e) 
0(e 3 ) 0(e) 0(1) 



(81) 



Again, we need to diagonalize the (2, 3) submatrix with the matrix W defined in Eq. (168)) 
with the angle defined in Eq. ()b9jl . The resulting matrix is H'" given in Eq. ()7Ujl . From 
this point on, we must treat the Sm^ > and Sm^ < cases separately since level crossing 
between a and Sm^ occurs for the Sm^ > case but not for the Sm^ < case. 

When 0777.3! > 0, we can use 1 — s' 12 = 0(e 4 ), and A 2 = ac 2 3 + 0(e 2 \5m\ l \) to approximate 



tan ! 



a sin 2#i3 



£777,01 — acos2#i 



+ 0(e 3 ) . 



(82) 



&3! — tx iuo ^-l/^3 

In this case, we expect = 0(1), = 0(1), and = 0(1). The A"s can also be expanded 
in 5m 21 /a = 0(e 2 ) and we find 



K = A'! 



(a + 5m 2 x ) — \J(a — 5m 2 31 ) 2 + 4a 5m 



2 „2 
31*13 



A3 — M 



+ 0V773-J + >/(a — 0V77IJ 2 + 4a Sm^s^ 



+ 0(<Kl - a)0777^ lS f 2 + 0(^^771^1) 



+ 6(a - 5mi x )5rnf 21 sl 2 + O^^m^l) 



(83) 

where 6 is the Heaviside step function. 

For the Sm^ < case, the denominator in the definition of tan 20 in Eq. (joTJj) is always 
negative and never crosses zero for any value of a. Therefore, the smallness of sin 2# 13 in the 
numerator is never cancelled by an equally small denominator. This allows us to expand 4> 
in #i3 and we find 

*^K+ < £3) = 0(£ >- < 84) 



which is actually valid for all values of a when 5m 31 < 0. Therefore, = 0(e) and 
= 0(1). The A"s are expanded as 



A" 



>+ ~ ac 2 13 + Sm 2 21 s 2 12 + 0(e 4 \5m 2 31 \) 
A" « 5m^ + as? 3 + 0(e 6 |5m^|) . 



(85) 



For both the > and Sm 3l < cases, both A 2 and A3 are of the same order as 



IJmjj-J. Therefore, 



H" 



K 

— aC 12 Ci3Si3S ( , 



ac 12 c 13 s 13 



^2 



A 3 



1 5m 



0{e 2 ) 0(e 3 - 4 ) 0(e 3 ) 
0(£ 3 < 4 ) 0(1) 
0(e 3 ) 0(1) 



(86) 

where the sizes of the (1, 2) and (2, 1) elements depend on whether = 0(e) (5m 31 < 0) or 
4> = 0(1) (5m 2 31 > 0). In either case, further diagonalization is not necessary. 

If we relax our accuracy requirement and allow for errors of O (e 2 \ Sm^ \ ) , then the A"s for 
both the Sm 31 > and 5m 2 31 < cases can be approximated by 



X'i « 



(a + Sm^) ± a/(o — Smii) 2 + 4a ^m^sf 



13 



(87) 



We will argue that this approximation is sufficient later. The eigenvalues are then: 



Ax « a; 



r 2 2 

dm 21 c 12 , 



A 2 « A'' « 



13 



+ (Sm^) ± i/(a — (JmU 2 + 4a tfmj^s 



2 

13 



(88) 



where the upper sign corresponds to the Sm^ > case, and the lower sign corresponds to 
the Sm^ < case. 



AS a/\5m 2 31 \ = 0(e~ l 



The relative sizes of the elements of H' in this case are 



H' 



ac 12 c 13 



ac 12 s 12 c\ 3 



aci 2 s 12 cl 3 as\ 2 c\ 3 



8ml x 



aci 2 ci 3 si 3 
asi 2 ci 3 si 3 



oci 2 c 13 si3 asi 2 Ci3Si3 as 



13 



+ 5m 31 



\5m 



0(e^) 0(e^) 0(1) 
Oie- 1 ) O(e^) 0(1) 
0(1) 0(1) 0(1) 



(89) 



As in all the previous cases, the (1,2) sub matrix is diagonalized by V defined as in Eq. (J44|) 
with tp defined as in ([45)1 . The partially diagonalized form is H" in Eq. ()46[). with A' 1; A 2 , 
and 9' 12 defined as in Eqs. fl4T|) and (f4*%)). The expansions of these quantities are: 



A; = 5m 2 21 c 2 12 + 0(e 5 \5m 2 31 \) = 0(e 2 \5m 2 31 \) , 

A 2 = ac 2 13 + Sm 2 21 s 2 12 + 0(e 5 \Sm 2 31 \) = O^-^m^l) 



(90) 



and 



0? 



12 



7T 

2 

6*12 ■ 



" #12 



2a 



sin(2# 



12 J 



Ofe £ 



7T 

2 



^iisin(2M+0(e 5 ) 
2a 



Since s' 12 = 0(1), and ac' 12 = 0(e 2 \5m 31 \), we find 









a; 



ac 12 C 13 s 13 

as' 12 ci 3 si3 



ac' 12 ci3Si 3 as' 12 ci 3 si3 as\ 3 



8m 2 31 



\8m 2 31 \ 



0{e 2 ) 


0(e 3 ) 







O^ 3 



0(0 0(1] 

0(1) orr 



(91) 



(92) 



We diagonalize the (2, 3) submatrix with the matrix W defined in Eq. (J68|) with the angle 
4> defined in Eq. (l69|) . The resulting matrix is H'" given in Eq. (J7Uj) . 

Recall that 20 is in the second quadrant if Sm 31 > (level crossing occurs), and in the 
fourth quadrant if 8m 31 < (no level crossing). The expansions of A 2 , A3, and </> in powers 
of 5m 3l /a = 0(e) differ accordingly. For the 5m 2 1 > case, we find 



\"_ 



6m 2 31 c 2 3 + 0(e 3 \6m 2 31 \) = 0(\6m 2 31 \) , 

a + 5m 2 31 s 2 13 + 5m 2 21 s 2 12 + 0(s 3 \5m 2 31 \) = 0(e- l \5m 2 31 \) 



(93) 



and 



2 



5m 2 31 



7 13 



0(e 3 ) = 0(1) , 



in which case both and are of order 1 and 



H" 



-ac' 12 Ci3Si3S^ 



-aC 12 Ci3Si3S (? i aC 12 Ci3Si3C<^ 











\5m 



For the Sm 31 < case, we have 



0{e 2 ) 0(e 3 ) 
0(e 3 ) 0(1) 
0(e 3 ) 



(94) 



0(e 3 ) 


Oie- 1 ) 



(95) 



6m 2 31 s 2 3 + Sm 2 21 s 2 12 + 0(e 3 \Sm 2 31 \) = O^l^il) 



A 3 



A'! = 5ml lC \z + 0{e"\5ml l \) = 0{\5m 2 x \) 



(96) 



and 



-9 



5m 



13 



31 Q 



13 



0(r 



0(e) 



(97) 



which is just Eq. ()84|) expanded in powers of 5m 31 /a. In this case, = 0(e) and = 0(1). 
Therefore, 



H" 



-ac' 12 ci3Si 3 S0 



ac 12 ci 3 si 3 c^ 



-aC 12 Ci3Si3S^, CLC 12 Cx 3 Si 3 C^ 



A 2 



A 3 



1 8m 



0(e 2 ) 0(e 4 ) 0(e 3 ) 
0(e 4 ) O^ 1 ) 
Ofe 3 ) 0(1) 



(98) 



In either case, the Hamiltonian has been approximately diagonalized. The eigenvalues for 



the Sm 31 > case are 



Ai 
A 2 

A 3 



a; 
a': 

A" 



5m 



2 2 
21 c 12 ) 

8m 2 31 c\ 3 , 



a + <5m 31 s 13 + Sm 21 s 12 



(99) 



while for the Sm 31 < case, they are 



Ai 
A 2 
A 3 



A'i 

K 
A" 



<5?7i 21 c 12 , 



5m 2 31 c 2 13 



21 c, 12 > 



(100) 



B. Effective Mixing Angles 

To summarize what we have learned above, when a/l^m^l = 0(e 2 ) or smaller, H' can 
be approximately diagonalized by a single (1, 2) rotation using the matrix 



V 



1 



(101) 



where 



with < <p < 
using the matrix 



|. When a/\5m 2 



311 



0(e) or larger, this must be followed by a (2, 3) rotation 



W 



where 



Ccj, 



cos i 



sin i 



tan 20 



1 


-s^ 



as' 12 sin 2^ i3 



(103) 



a sin 2$ 



13 



+ as\ 3 



A' 



(Jm^ — a cos 2#i3 



(104) 



with s' 12 and A' 2 defined in Eqs. (g7D and (@HD- If Sm 2 31 > 0, then < < § . If Jm^ < 0, 
then -f < < 0. 

Note that the first case is encompassed in the second, since in the first case the sec- 
ond rotation angle becomes negligibly small. Therefore, the matrix which approximately 
diagonalizes the effective Hamiltonian is 



U' 



UQVW 

1 

C23 S23 

-S23 c 23 

1 
C23 S23 

— S23 C23 

1 
C23 S23 
— S23 C23 



C13 s 13 e~ ld 

1 

si 3 e t5 C13 

C13 s 13 e~ lS 

1 

-s l3 e lS ci 3 

C13 s 13 e- 4<5 

1 

4 <5 



-sise* 



C13 



C12 S12 

S12 C12 

1 

C12 S12 

-S12 C12 

1 



C19 s^2 



-12 



'12 u 12 



1 



1 
1 

e iS 



s v 



s v 



1 

1 

1 

e iS 



1 

1 
1 

e iS 

1 


—s^ 

c 13 s 12 s + s 13 c <f> 



w 



C13C12 c 13 s 12 c 4> ~ s 13 s 4> 

~ c 23 s 'l2 — • s 23 s 13 c/ 12 e * <5 c 23 c 'l2 C </> — s 23 ( s 13 s/ 12 C ^ + c 13 S (/>) el5 c 23 c 'l2 S — S 23( S 13 S 12 S — C \3 C ^)^ 
S23S12 - C23Sl3C / 12 e J<5 -S 2 3ci 2 C0 - C 23 (Si3si 2 C^ + C^S^e 15 -323^25^ - C 2 3(Sl3s'l2 s <?i - C 13 C^)e l& 

(105) 



We would like to identify this matrix with 



U 



1 

C23 S23 
— S23 C23 



C13 



si 3 e 



15 



srse" 1 

1 

C13 



C12 S12 
-S12 C12 
1 



a 




c 'l2 






c 13 


(1 - s' 12 )s^ 


si3c' 12 (l - s' 12 )s^ 


\ c 13 / 


C '\_2 S 4> 


c' 12 S0/c' 13 


lr\m 2 1 

|0m TI | 






0(e 3 ) 


± 


O(l) 


O(l) 


0( £ 4 ) 


O(l) 


0( £ 4 ) 


0(s 5 ) 


0( £ 5 ) 


0(e 4 ) 


0(e 4 ) 


0(e 2 ) 


± 


O(l) 


O(l) 


0( £ 3 ) 


O(l) 


0(e 3 ) 


0(£ 4 ) 


0( £ 4 ) 


0(£ 3 ) 


0(e 3 ) 


0(e) 


± 


0(e) 


0(e 2 ) 


0(e 2 ) 


O(l) 


0(e 4 ) 


0(e«) 


0( £ 4 ) 


0( £ 3 ) 


0( £ 3 ) 


O(l) 


+ 


0(e 2 ) 


0( £ 4 ) 
0( £ 4 ) 


0(e) 
O(l) 


O(l) 
O(l) 


0(£ 5 ) 
0(e 4 ) 


0(£ 8 ) 
0(£ 7 ) 


0( £ 4 ) 


0(e 3 ) 
0(e 2 ) 


0(e s ) 
0( £ 2 ) 


OCe- 1 ) 


+ 


0( £ 3 ) 

0( e 3 ) 


0(e 6 ) 
0( £ 6 ) 


0(e) 
O(l) 


O(l) 

o( £ 2 ) 


0(e 7 ) 
0(e 4 ) 


0( £ 8 ) 


0(e 5 ) 


0( £ 4 ) 
0( £ 3 ) 


0(e 4 ) 
0(e) 



TABLE II: The sizes of the factors (1 
and cf 12 s^/cf 13 . 



C12C13 



'12 



)^/ c 'l3> s 13ci2(! - s 'l2) s 0/ c 'l3' ^12 C 1 - C 'l3/ C 13) 5 c' 12 S0, 



S12C13 



„i<5 



-S12C23 - Ci 2 Si 3 S 2 3e C12C23 - Si 2 Si 3 S 2 3e C13S23 



(106) 



S12S23 ~~ Ci2Si3C23e — C12S23 — Sl2 s 13 c 23 e C13C23 

up to phases that can be absorbed into the Majorana phases of the neutrinos and redefinitions 
of the charged lepton fields. Comparing the (1, 3) elements of the matrices, we can make 
the identification 



Sl3 



c 13 s 12 s 4> + s 13 c 

(Cl 3 S0 + - (1 - S' 12 )S0C1 3 



'13 



13 



(1 - s'^s^cn j 



(107) 



13 



where we have defined 



s' 13 = sin 0' 13 , c' 3 = cos ^ 3 , = 9 13 + 



'13 ) "-13 — ^di7 13 , (7 13 — (713 I" ^ • (108) 

Note that the factor (1 — s'^s^/ c! 13 is of order e 3 or smaller regardless of the value of a, as 
shown in Table HU Therefore, 



which implies 



sin ^13 = sin [9' 13 + 0(e c 



9 13 = 9' 13 + 0(e 3 ) . 



(109) 



(110) 



Next, looking at the (1, 1) and (1,2) elements, we find 

tan# 12 = ^ 12 f 13 
C12C13 



c 13 s 'l2 C </> ~ s 13 s (j> 
C 13C 12 



(C13C0 — Si 3 S^)s' 12 — S 13 (l — S 12 )S^ 



Cl3Ci 2 



c 13 



taufl, ' ' ^12(1-^)^ 



12 co S 2 9' 12 1 



111^ 



From Table ITT1 we find that the factor s 13 c' 12 (l — s' l2 )s^ ) /c' l3 is of order e A or smaller for all 
a. Therefore, 



tan# 12 = \^jtan[9' 12 + 0(e A )} . (112) 

Since 9' 13 = 9i 3 + (j), we can expect the ratio c' 13 /ci 3 to be roughly equal to one when <p is 
small, and consequently, 9±2 ~ 9' 12 . Indeed, if 5m 31 > with a/l^m^l < 0(e), or 5m| x < 
with any a, then < 0(e), and we find 



1- — = (I - c^) + s^tan9 13 < 0(e 2 ) . (113) 

Cl3 

In these cases, we can treat (1 — c' 13 /ci 3 ) as a small quantity and expand 

9 l2 = 9' 12 + s' l2 c' l2 (l - &\ + ■ - - = 9' 12 + 0(e A ) . (114) 

The 5m 31 > case with a/\Sm 31 \ = 0(1) or a/l^mU = 0(e~ 1 ) must be considered sepa- 
rately. First, taking the reciprocal of both sides of Eq. ()112|) . we obtain 

cot 0i2 = ^ cot 9' , (115) 

Cl3 

where we have dropped the shift in 9' 12 on the right hand side which is of order e 7,8 in these 
particular cases, (cf. Table ITT1) Recall that when a/l^mg-J = 0(e°'~ 1 ), we have 

^2 = f-^sin(20 12 ) + O(^), (116) 



while 



Therefore, from Eq. ()115|) we find 



^=0(e°'~ 2 ). (117) 

c 13 



( ? - 6i2 ) = ~t~ tan ( tt ~ &i2 



2 c' 13 \ 2 



^tan ( 6 -^sin(29 12 ) + 0(e 
c 13 V 2a 



tan f^^iisin(2M + 0(e 4 ' 3 ) 



■ 



from which we can conclude 



7T C13 / , 



Next, using the relation 

Sl3s'i2 S <^ _ C 13 c <£ = -(C13C^, - - Si 3 (l - S' 12 )S0 = -c' 13 + 0(e 4 ) , 

we simplify the (2, 3) and (3, 3) elements of U' as 

C23C' 1 2 S </' ~ •5 2 3(S13S / 12 S <} « ) - C^C^e* 5 = C 2 3c' 12 S <j) + S 2 3C 13 e %S + 0(e 4 ) , 

-S23c'i2 s J> ~ c 2 3(si3s' 12 s (f) - c^c^e 15 = -s 23 ci 2 s + c 23 c' 13 e lS + 0(£ 4 ) 
Then, using the fact that c' 12 s^ = 0(e 2 ) or smaller, we find 
\c23c'i2 s <t> - ^(s^s'^s^ - ci 3 c<p)e zS \ 



2 + 2s 23 c 23 c / 13 c / 12 s </j cos 6 + c 2 23 c' 12 2 sl + 0(e 4 ) 



S 23 C 13 



s 23 c 'l3 + 2S23C23C / 13 C / 12 S0 COS 5 + 0(e 4: ) 
S23c'i 3 + C2 3 c' 12 S0 COS 5 + 0(s^) 
s 23 c 'l 2 s <j> — c 23( s 13 s 'l2 s <^ — c 13 c cf>)G lS \ 

c 23 c 'i3 2 - 2 s 23 c 2 3C , 13 c' l2 s lj) cos 5 + s 2 23 c' l2 ! sj + 0(e 4 

C 23 C 'l3 2 _ 2s 2 3C23C].3Ci2 s ^ COS ^ + O^ 4 ) 
C23Ci 3 - 523^50 COS 5 + 0(e 4 ) 



Therefore, we can make the identification 

Cl 3 S 2 3 



tan 9 



2:-! 



C13C23 

— r 1 + °( e 

C23Ci 3 - S 23 C 12 S0COSO 
t 23 +(^t) COS 5 

l-t2 3 (^)cOS, 



tan 



e 23+ (^t) cos 5 

V C 13 / 



0(e 



4n 



and we obtain 



c 13 



6 23 = 6 23 +[ C -^)cos5 + 0(e*). (124) 



The factor d 12 s^/ c' 13 is of order e 3 or smaller if 8m\i > with a/IJmjjJ < 0(e), or 5m 31 < 
with any a. In those cases, we have 

e 23 = 6 23 + 0(e 3 ) . (125) 
For the case of 5m 31 > with a/IJmjj-J > 0(1), we can expand c' l2 and approximate 

23 = #23 + 4^ sm(2^ 12 ) cos 5 + 0(e 4 ) . (126) 



c' 13 \ 2a 

Finally, we calculate the CP violating phase. The Jarskog invariant of U' is 

J' = (c 13 s' 12 s^ + s 13 c^)(cvis' 12 c^ - s 13 s^) (c 13 c' 12 )s 23 c 23 sin 5 

= Si3(Ci3Si2)(Ci3Ci2)s 23 C23Sin(5 

^ ^2 ^ 

= (si3C 13 Si2Ci2)s23C 2 3Sin5 . (127) 

On the other hand, the Jarskog invariant of U is 

■J = s 13 c 13 si 2 c 12 s 23 c 23 sin 5 . (128) 

Comparison with J' shows that 

sin(2^2 3 ) sin 8 = sin(2# 23 ) sin 5 , (129) 

which is actually an exact relation as discussed in Ref. Since #23 = 6*23 + 0(e 3 ) when 
Sm 31 > with a/l^mtjjj < 0(e), or 8m 31 < with any a, for these cases we have 

S = 5 + 0(e 3 ). (130) 

For the case 8m% x > with a/l^m^] > 0(1), we can use Eq. (|12fij) to obtain 

5 = 5-^( 5 -^) ^4 sin 5 + 0(, 6 ) . (131) 
c' 13 \ a J tan(2fl 23 ) v ' K ' 

C. Summary of Neutrino Results and Sample Calculation 



Let us summarize the results of the two previous subsections. 



Approximate values of the effective mixing angles in matter can be obtained from the 
relations 



013 ~ 013 > 

tan « — tan9' 12 , 

Cl3 

023 « 23 +(^)cOs5, 

sin(20 2 3)sin5 = sin(20 23 ) sin 5 , (132) 



where 



012 = 012 + V 2 , 

01 3 = 013 + 0, (133) 

and the angles ip and were defined in Eqs. (jlKjl and (jflEJ), respectively. The a-dependence 
of ip and for the sample case of tan 2 i2 = 0.4, sin 2 (20i 3 ) = 0.16, Sm^ = 8.2 x 10~ 5 eV 2 , 
and l^m^l = 2.5 x 10~ 3 eV 2 is shown in Fig^with gray solid lines. In the figure, the angles 
are in units of ir, and they are plotted against the variable a defined as: 

a " log ^ wk\ ' wk\ = £ ~ a • (134) 

a = corresponds to a = ISm^], and a = — 2 corresponds to a = ^m^. 

When either Sm^ > (normal hierarchy) with a < — 1, which corresponds to a/l^m^l < 
0(e), or 5m| x < (inverted hierarchy) with any a, the angle is small, and Eq. ()132|) reduces 
to 

013 ~ 013 = ^13 + ) 
012 ~ 012 = 012 + V 3 ; 
023 ~ 023 , 

5^5. (135) 

For the 5m 2 ;! > case (normal hierarchy) with a > 0, which corresponds to a/l^m 2 ^ > 0(1), 
the angles can be approximated as 




FIG. 1: The exact (gray solid line) and approximate (black dashed line) values of tp and eft plotted 
against a = \og l i e {a/\5m\ l \). The parameter choice was tan 2 6\2 = 0.4, sin 2 (2#i3) = 0.16, Sm^i = 
8.2 x 10~ 5 eV 2 and |5m^| = 2.5 x 10~ 3 eV 2 . 




FIG. 2: Aip = 93 a pp r ox - ^cxact and A0 = ^ apP rox - '/'exact plotted against a = log 1/e (a/|(5m| 1 |) for 
the same parameter choice as Fig. ^ 



23 « 02 3 + ^( 5 -^) sin(20 12 )cos5 
c' 13 V 2a 7 



c' 13 \ a J tan(20 23 ) V ; 

To make use of these expressions, we must first calculate if and 0, and then 9' 12 = 12 + 93 
and 6^ 3 = ^13 + 0. Simple approximations to and are provided by 



a sin 26*19 a sin 20i 3 

tan2<^ « _ **— , tan 20 « _ . 137 

" m 2i — a cos 2^12 tojj — a cos 2t7i3 



The approximate values of (p and obtained from these expressions are also shown in Fig. ^ 
with black dashed lines. As is clear from the figure, the graphs of the exact and approximate 
values are virtually indistinguishable at this scale. In Fig. El we plot the differences between 
the approximate and exact values of tp and 0: 

A(f = ^approx — V^exact j A0 = 0approx — 0cxact • (138) 

The maximum deviation from the exact values occur at the level- crossing points a = —2 
and a = (there is no level- crossing at a = when 8m 31 < 0) but even then, it is a mere 
fraction of a percent of tt. Using Eq. (|137jl . we can also obtain simple approximate formulae 
for 9' 12 = 9 12 + tp and 9[ 3 = 13 + 0: 

tan2e ; 2 „ , ta„2#; 3 « Jjpk^hi . (139 ) 

om 2l cos 26 , 12 — a " m 3i cos 2#i3 — a 

These allow us to calculate 9[ 2 and 9' 13 directly without going through <p and 0. Eqs. (|137|) 
and (|139j) provide a quick and easy way to obtain the input angles necessary to utilize 
Eqs. (HS2D, (gnilD , and (fTT^jt . 

To demonstrate the accuracy of these approximations, we present a sample calculation 
using the following parameter choice: 

5m 2 21 = 8.2 x 10" 5 eV 2 , 



\5m 2 31 \ = 2.5 x 10" d eV 2 , 

tan 2 = 0.4 , 

sin 2 (20 13 ) = 0.16 , 

023 = 0.2 7T , 

S = 0.25 tt . (140) 

The values of 5m 2l , \5m 31 \, and tan 2 12 are the experimental central values. The value of 
sin 2 (20 13 ) is taken to be the 90% upper limit corresponding to our choice of l^mU so that 



the 0(6* 13 ) terms that we neglect are maximized. The value of 6*23 is also chosen to be one of 
the 90% confidence limits since if we set 623 to the experimentally preferred central value of 
0.257T, then the shift of 5 away from 5 would be suppressed. Similarly, we chose 5 = 0.257T 
so that both 823 and S will be shifted from their vacuum values. 

We first consider the normal hierarchy case (Sm^ > 0). In Fig. EH, we plot the exact 
values of #12, # 13 , 6*23, and 5 calculated numerically with gray solid lines, together with the 
approximate values obtained from Eq. (|132|) . using Eqs. (|137|) and (|139|) to calculate the 
input angles, with dashed black lines. Fig. Eh shows the errors: 

A9. L j = 9ij t approx exact , AS = 5 a pp rox ^cxact • (-^-0 

A0 12 and A613 are indicated with solid gray lines, while A# 2 3 and A5 are indicated with 
dashed black lines. The errors are never larger than a fraction of a percent of n, and 
comparison with Fig. E] makes it apparent that the majority of it was inherited from having 
calculated ip and using Eq. (jl37J) . In Figs. EH and EH, we plot the exact values against 
the approximate values obtained using Eqs. (j!35|) and ()136|) . respectively, together with the 
errors in Figs. Eb and Eh- Clearly, the approximations using Eqs. (jl35)l and (|136j) are good 
in their respective ranges of applicability. 

For the inverted hierarchy case {Sm^ < 0), we only need to consider Eq. ()135|) . In 
Fig. EH, we show the comparison between the numerically calculated exact values and the 
approximate values obtained from Eq. (jl35|) . The errors are shown in Fig. Eb- 

The approximate values for the mass-squared eigenvalues are given by 

Ai w A'_ , 

A 2 w a': , 

A 3 « A'; , (142) 

for Sm^ > (normal hierarchy), and by 

\i « A'_ , 

A2 ~ A'| , 

A 3 « A'! , (143) 

for Sm^ < (inverted hierarchy), where \'± and A± are defined in Eqs, ()49j) and (|72jl . 
respectively. The accuracy of this approximation is illustrated in Fig using the parameter 



Normal Hierarchy Normal Hierarchy 




-3-2-1012 -3-2-1012 
a a 



FIG. 3: (a) The exact values of 612, 613, 623, and 5 (solid gray lines) plotted against their 
approximate values (black dashed lines) obtained using Eq. ()132j) . with Eqs. (|137j) and ()139|) . as 
functions of a = log 1/£ (a/\5mj 1 1). (b) The differences A9 12 = 6>i2, apP rox - #i2,exact and A6»i 3 = 
^i3,a PP rox - ^i3,cxact (solid gray lines), and the differences A6> 23 = #23,a PP rox - #23,cxact and A5 = 
<5approx — 5 C xact (black dashed lines) of this approximation plotted against a. 
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FIG. 4: (a) The exact values of 9\2, 613, 623, and 5 (solid gray lines) plotted against their approx- 
imate values (black dashed lines) obtained using Eq. (|135|) as functions of a = log 1 / £ (a/|(5m| 1 |). 
(b) A$i2 and A613 (solid gray lines), and A#23 and A5 (black dashed lines) of this approximation 
plotted as functions of a. This approximation is applicable when a < — 1. 
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FIG. 5: (a) The exact values of #12, #13, 023, and 5 (solid gray lines) plotted against their approx- 
imate values (black dashed lines) obtained using Eq. (|136l) as functions of a = logy e (a/\Sm^i\). 
(b) A#i2 and A#i3 (solid gray lines), and A#23 and A 6 (black dashed lines) of this approximation 
plotted as functions of a. This approximation is applicable when a > 0. 
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FIG. 6: (a) The exact values of 9\2, 613, 623, and 5 (solid gray lines) plotted against their approx- 
imate values (black dashed lines) obtained using Eq. ()135|) as functions of a = log 1 / E (a/|5m3 1 |) f° r 
the inverted hierarchy case (Sm^ < 0). (b) A#i2 and A#i3 (solid gray lines), and A#23 and A5 
(black dashed lines) of this approximation plotted as functions of a. 
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FIG. 7: The exact and approximate values of log 1 / e (5A2i/|^| 1 1) and log 1 / e (|5A3i|/|(5m3 1 |) for the 
parameter set of Eq. (|140|) plotted against a = log 1 ^ £ (a/\5m'^ 1 \). The exact values are in gray solid 
lines, whereas the approximate values are in black dashed (<5A2i) and black dot-dashed ( | <5A3i [ ) 
lines. 
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FIG. 8: The rescaled errors AA21 (dashed) and AA31 (dot-dashed), as defined in Eq. ()144)) . for the 
approximation of Fig. [7J 
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FIG. 9: Comparison of exact and approximate values using Eq. (|146f) for the normal hierarchy 
case. The approximation is applicable when a < — 1. 
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FIG. 10: Comparison of exact and approximate values using Eq. 1)147(1 for the normal hierarchy 
case. The approximation is applicable when a > — 1. 
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FIG. 11: Comparison of exact and approximate values using Eq. ()146() for the inverted hierarchy 
case. 



values of Eq. (|14(J|) . where the exact numerically calculated values of <5A 2 i = A 2 — Ai and 
I^Asi | = | A 3 — Ai| are plotted against those obtained from the above approximate expressions. 
The vertical axis is plotted using the same log-scale as the horizontal axis where l^m^l 
corresponds to and Sm 2 -^ corresponds to —2. Since a log-scale plot does not reflect the 
absolute accuracy of each SXij, in Fig. |S] we plot the difference between the exact and 
approximate values of SXtj normalized to £A miniCxact : 

A \ oAji approx oXn exact / • o o\ I a a a\ 

AAa = — , (z = 2, 3) , (144) 

A m i nGxac t 

where 

5A m i njGxact = min(5A2i,cxact, | <5Av3i ;exact | , |<5A3 2iex act|) • (145) 

This tells us how large the errors are compared to 5A m j n , which is typically used to expand 
the oscillation probabilities in. As is evident from the figures, the approximation is excellent. 

When either Sm^ > (normal hierarchy) with a < —1, i.e. a/lSm 2 ^ < 0(e), or 
Sm^i < (inverted hierarchy) with any a, the A's can be further approximated by 

(a + Sm 2 ,^ — a/(o — Sm^) 2 + 4a 6m\ x s\ 2 



Ai 

A 2 
A 3 



(a + Smli) + a/(o — Sm^) 2 + 4a Sm^s^ 



5m 2 31 . (146) 
For the Sm^ > (normal hierarchy) case with a > —1, i.e. a/lSm 2 ^ > 0(e), we can use 



Ai « Sm 2 21 c 2 12 , 
(a 

A 2 t 



a + dm^) — y(a — Sm^) 2 + 4a dm^s^ 



■ 



Ag _ (a + SmlJ + y/(a - (5m^) 2 + AaSm 2 ^ 

This second approximation introduces errors of O^l^m^l) in the A's. However, since 
SX21 > 0(£r|(57Ti§ 1 |) in this range, an error of this size is tolerable. The accuracy of these 
approximations is illustrated in Figs. El and El for the 5m^ 1 > case, and Fig. ^2 f° r the 
8m\ x < case. Though the accuracy is not as good as when Eqs. ()142|) and ()143|) . with 
Eqs. (jinj) and f72*)l . are used, it is sufficient for most practical purposes. 



D. Anti-Neutrino Case 



Matter effects for anti-neutrinos can be treated in an analogous fashion. We will therefore 
omit the details and give only an outline of the derivation and results. 
The effective Hamiltonian for anti-neutrinos in matter is 



(148) 





Ai 





















—a 








H = U* 





A 2 





U T = u* 





8m 2 21 





u T + 



















As 
























We denote the effective mass-squared eigenvalues as Aj 



1, 2, 3), and the diagonalization 



matrix as U to distinguish them from those for the neutrinos. (Note the mirror image of 
the tilde on top of U.) This matrix can be partially diagonalized as 



H' 



QU HU*Q* 




Sm 2 21 
8m 2 31 



Q 



y. L 



5m 2 



+ t/ J 



21 








5m 2 31 



Q*~aQ 



-a 

U* \ Q* 


u el u* el u el u* e2 u el u: 3 

u e2 u* el U e2 U* e2 U e2 U* e3 

u e3 u* el u e3 u* e2 u e3 u: 3 



Q* 





Sm^ 








5m 2 31 



2 ^2 



C 12 C 13 Ci 2 S 12 cf 3 C12C13S13 
C\ 2 S 12 C\ 3 S 2 l2 C 2 13 S12C13S13 



C12C13S13 S12C13S13 



'13 



-acToC' 



12^13 



-ac 12 s 12 cl 3 



-aci 2 c 13 si 3 



-aci 2 si 2 c{ 3 -as{ 2 c{ 3 + 5m 2l -asi 2 c 13 s n 

,2 



5m 2 31 



(149) 



-CZC12C13S13 — asi2Ci3Si3 — as 13 

The only difference from the effective Hamiltonian H' for the neutrinos, Eq. (|38|) . is in the 
sign in front of the matter-effect term a. 

Following the general procedure we employed for the neutrinos, we begin by diagonalizing 
the 1-2 submatrix of H'. Let 



V 



1 



(150) 



where 



ac\n sin 2#i 2 ( 

= cos <p , s v = sm<p, tan2y? = - 2 — — , - 

om^i + acj 3 COS 26 , i 2 V 



:~<^<o) . 



Using V, we find 



/J" = V^H'V = 



o 

— ac' 12 Ci3Si3 





A 2 

— as' 12 Ci3Si3 



-ac' 12 ci 3 si 3 

— aS 12 Ci3Si3 
-« S 13 



if 3 + 0V4 



tan2? 19 = fill^ 

om 21 cos 2^12 + ac( 3 



where 

c' 12 = cos 9' 12 , s' 12 = sin 9' 12 , 9' 12 = 9 12 + (p , 
with 

(grafi ~ acj 3 ) ± V^ro^ + ac 2 3 ) 2 - 4ac 2 3 5m 2 21 s{ 2 



and A' 



>i — A_, A 2 — A + 



Next, unlike the neutrino case, 



we 



diagonalize the 1-3 submatrix of H" . Let 



W = 



Ccf, S0 

1 

-s^ 



where 



The angle is in the fourth 
Using W, we find 



(151) 



(152) 



(153) 



(154) 



(155) 



7_7 - oc' 12 sin 2# 13 

c«* = cos0, SA = sm(p, tan20 = --— 5 =- . (156) 

5m 2 31 - as 2 l3 - X[ 

l the fourth quadrant when 8m 3l > 0, and the first quadrant when 5m 31 < 0. 
d 



H'" = W^H"W = 



as/ 12 C 13 s 13 s 



as 12 c 13 s 13 s ^ 

A' 2 



-as 12 i 



— as^ 2 Ci3Sx3C0 





: C 13 S 13C0 

// 

'3 



A" 



(157) 



where 



A" — A" , A 3 ' — A' + , 
A/ = A" , A3 = \"_ , 



if Sm^ > , 

if 8m 2 3l < . 



(158) 



with 



[{5m 2 , - as 2 ,) + \\ ± v / [(5m| 1 -a S ? 3 )-A' 1 ] 2 + 4a2c' 1 2 2 c? 3S 



13 



(159) 



Evaluation of the off-diagonal terms of H'" reveals that it is approximately diagonalized, 
and the diagonalization matrix is given approximately by 



U'* = U*Q*VW = U*VQ*W 



(160) 



or taking the complex conjugate, 



U' 



UVQW 

1 

C23 S23 

—s 23 c 23 

1 

C23 S23 

— S23 C23 



C13 s 13 e- i5 

1 

-s l3 e lS C13 

C13 s 13 e~ lS 

1 

45 



-8i 3 e' 



Cl3 



C12 S12 
-812 C12 
1 





12 



12 

^12 ^12 



C 13 C 12 C <^ ~~ s 13 s ^ 



1 



c 13<Sl2 







1 

1 
1 

e iS 



QW 



" c 23 s 'l2 C — s 23( c 13 s </> + s 13 c 'l2 c <?0 



■=1 „iS 



C 23 c 12 ~ s 23 s 13 s 12 e 



c<f, 
1 
-s<j, 

■513^0 + Ci3C^ 2 S^ 
C 23S' 12 S0 + S23(Ci3C^ — Si3C / 12 S (j / ) )e J ' 5 



S23s'i2 c </> - c 2 3(ci3S0 + siac^c^e* 5 -s 2 3c' 12 - c 23 s 13 s' 12 e ld s 23 s' 12 s <t) + C2 3 (ci 3 C^ - s 13 c' 12 s 



iS 



48 



(161) 



Identification of this matrix with 



U 



1 

c 2 3 S23 

U — S23 C23 

C12C13 



C13 s 13 e~ i5 
1 
-sise" ci 3 

S12C13 



C12 S12 
-S12 C12 
1 



"S12C23 — c 12"Sl3 s 23 e 



S 12S23 — Ci2Si3C23e 



C12C23 — Sl2"Sl3 s 23 e 



i<5 



"C12S23 — s 12 s 13 c 23 e 



i<5 



si3e" 



C13S23 



C13C23 



(162) 



(up to phases that can be absorbed into redefinitions of the charged lepton fields and the 
Majorana phases of the neutrinos) yields 



9 



13 



tan 6*12 

#23 



13 ) 
Cl3 



tan B\ 



12 5 



13 



#23 — 



/ s' 12 s<t> 



13 



cos 5 



sin(20 23 ) sin 5 « sin(20 23 ) sin 5 , (163) 

where 



= sin 9' 13 , c' 13 = cos 9[ 3 , 6' 13 = #i 3 + <fi . (164) 



'13 



For the 8m 31 > case (normal hierarchy) with any a, or the 8m 31 < case (inverted 
hierarchy) with a/l^m^] < 0(e), the angle is small and the above relations simplify to 

#i3 « e' 13 = e 13 + , 

#12 ~ #12 — #12 + <f , 



#23 ~ #23 



5^5. (165) 

For the < case (inverted hierarchy) with a/IJm^l > 0(1), we can approximate 



'23 



#23-^f^)sin(2Mcos5, 



c' 13 V 2a 



- j- , (8m 2 21 \ sm(26»i 2 ) . . 

5^5 + -^ — v y sm 5 . (166) 

c' 13 V a y tan(20 23 ) 

The behavior of the angles tp and 0, defined in Eqs. ()151j) and ()156|) . are plotted in 
Fig. ^1 as functions of a = log 1 / £ (a/\Sm 3l \) with gray solid lines for the parameter choice of 
Eq. Ijl4()j) . Approximate values can be obtained from 

a sin 2#i2 - asin20i 3 , . 

tan2yj^-— - -5 — — , tan20^-— - — — , (167) 

077% + Ct COS 2#i2 " m 31 + a cos ^#13 

which are also shown in Fig.^Jwith black dashed lines. The accuracy of this approximation 
is shown in Fig. (|13jl . The approximate values of 6[ 2 = #12 + <p and #'i 3 = #13 + can be 
calculated from the expressions 

+ 5m^sin2#i 2 =•, om^ sin 20 i3 

tan 20, 9 = - — ^—^ , tan 20, , = - — ^— — . 168 

om^i cos 20i2 + « om^ cos 20i3 + a 

The accuracy of Eqs. ()163J) . ()165|) . and ()166|) . using Eqs. (|167|) and ()168|) as input, is illus- 
trated in Figures IT41 through 1171 for the parameter set of Eq. (|140J1 . 

The approximate values for the effective mass-squared eigenvalues are given by 



Ai » \'L , 




FIG. 12: The exact (gray solid line) and approximate (black dashed line) values of ip and </> 
plotted against a = ^ogu e (a/\Sm^i\). The parameter choice was tan 2 6»i2 = 0.4, sin 2 (26»i 3 ) = 0.16, 
5m 2 2l = 8.2 x 10~ 5 eV 2 and \8ml x \ = 2.5 x 10~ 3 eV 2 . 



. 0005 


-0 . 0005 
K -0.001 
-0 . 0015 
-0 . 002 



Normal Hierarchy 



-3 -2 



-1 
a 



















Aip 
























A&/ 

































. 006r 

. 004 - 
.002 - 
: 

-0 .002 - 
-0 . 004 - 



Inverted Hierarchy 

































-J 






















A(p 





















-3 -2 



-1 
a 



FIG. 13: Aip = ^ a pp r ox - ^exact and Acj) = ^ a pp r ox - (/'exact plotted against a = log 1 / £ (a/|(5m| 1 |) for 
the same parameter choice as Fig. 1121 
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FIG. 14: (a) The exact values of 6u, #13, 623, and <5 (gray solid lines) plotted against their 
approximate values (black dashed lines) obtained using Eq. ()163|) . with Eqs. (|167|) and ()168|) . as 



functions of a = log 1 / e (a/\Sm^ 1 \). (b) The differences A812 = 6 , i2, a pprox — 



12, exact 



and AO 



13 



#i3,a PP rox - #i3,cxact (solid gray lines), and the differences Af9 23 = 6*23,a PP rox - #23,cxact and AS 



^approx — ^cxact (black dashed lines) of this approximation plotted against a. 
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FIG. 15: (a) The exact values of 6*12, #13, #23, and S (gray solid lines) plotted against their approx- 
imate values (black dashed lines) obtained using Eq. (|165|) as functions of a = log 1/ / £ (a/|5m| 1 1). 
(b) AQ\2 and A#i3 (gray solid lines), and A#23 and AS (black dashed lines) of this approximation 
plotted as functions of a. This approximation is applicable when a < — 1. 
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FIG. 16: (a) The exact values of 612, 9 13, #23 5 and 5 (gray solid lines) plotted against their approx- 
imate values (black dashed lines) obtained using Eq. (|166|) as functions of a = log 1 / £ (a/|5m3 1 1). 
(b) A#i2 and A^i3 (gray solid lines), and A#23 and A5 (black dashed lines) of this approximation 
plotted as functions of a. This approximation is applicable when a > 0. 
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FIG. 17: (a) The exact values of 6*12, #13, #23, and 5 (gray solid lines) plotted against their approx- 
imate values (black dashed lines) obtained using Eq. 1)165(1 as functions of a = log 1/ / e (a/|<5m3 1 1) 
for the normal hierarchy case (Sm^ > 0). (b) A^i2 and A^i3 (gray solid lines), and A023 and A5 
(black dashed lines) of this approximation plotted as functions of a. 



A 3 « K > (169) 
for the 5?n 2 ^ > case (normal hierarchy), and 

Ai « A'j , 
A2 ~ A' + , 

A 3 « A" , (170) 

for the Sm^ < case (inverted hierarchy), where A' ± and A± are defined in Eqs. ()154|) and 
f)159p . respectively. The accuracy of this approximation is illustrated in Fig [TH] using the 
parameter values of Eq. (jl4()|) . where the exact numerically calculated values of 8X21 = A2— Ai 
and \5\31\ = I A3 — Ai| are plotted against those obtained from the above approximate 
expressions. In Fig. EH we plot the difference between the exact and approximate values of 
SXij normalized to 5A min>exact : 

A T oXn approx Aji exact / . n \ ; 17 ,-, 

AAa = -= , (1 = 2, 3) , (171) 

min, exact 

where 

<5A m in jC xact = rnin(5A21,exact) |<5A3i i0x act|, 1 5 A32, ex act | ) • (1?2) 

When either 8rn^ 1 < (inverted hierarchy) with a < —1, i.e. a/IJm^l < 0(e), or 
Sm^ > (normal hierarchy) with any a, the A's can be further approximated by 

(5rri2i — a) — a/ (Sm^ + a) 2 — 4a dm^sf 



Ai 



'12 



2 

(5m 2 ,! — a) + \J (6mYi + a) 2 — 4a dm^s^ 



A 3 « . (173) 

For the c^m^ < (inverted hierarchy) case with a > — 1, i.e. a/l^mU > 0(e), we can use 



(Sm^i — a) + a/ (5^3i + a) 2 — 4a 5m 



2 „2 
31 6 13 



X 2 « bm\ x c\ 2 , 

_ (5m^ - a) - a/ (5m 2 ! + a) 2 - 4a fon^s^ 

The accuracy of these approximations is illustrated in Figs. [20 and |^ for the dm^ < case, 
and Fig. 1221 for the 8m\i > case. 
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FIG. 18: The exact and approximate values of log 1 / e (5A2i/|^3i I) and log 1 / e (|5A3i|/|(5m| 1 |) for the 
parameter set of Eq. (|14U|) plotted against a = log 1 / e (a/|5m3 1 |). The exact values are in gray solid 
lines, whereas the approximate values are in black dashed (SX21) and black dot-dashed (|#A3i|) 
lines. 
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FIG. 19: The rescaled errors AA21 (dashed) and AA31 (dot-dashed), as defined in Eq. (|171|) . for 
the approximation of Fig. 1181 




FIG. 20: Comparison of exact and approximate values using Eq. 1)173(1 for the inverted hierarchy 
case. The approximation is applicable when a < — 1. 
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FIG. 21: Comparison of exact and approximate values using Eq. (|174|) for the inverted hierarchy 
case. The approximation is applicable when a > — 1. 
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FIG. 22: Comparison of exact and approximate values using Eq. (|173() for the normal hierarchy 
case. 



IV. SAMPLE CALCULATION OF OSCILLATION PROBABILITIES 





L (km) 


P (g/cm 3 ) 


E (GeV) 


Q = lo gi/ e ( a /l 5m iil) 


|Asi| 


Reference 


T2K (JPARC -> Supcr-K) 


295 


2.6 


0.25 ~ 2 


-2.3 1.1 


(0.3 ~ 2.4)?r 


110] 


JPARC -> Korea 


1000 


2.7 


1 ~ 6 


-1.5 0.4 


(0.3 ~ 2.0)tt 




BNL — > Home Stake 


2540 


3.4 


2 ~ 10 


-0.9 ~ 


(0.5 ~ 2.6)tt 


Oil 



TABLE III: The three cases for which we calculate the probabilities for the processes Vn — > Va 
and i>u — ► z/ e . The ranges of a and IA31I were calculated assuming Sm^x = 8.2 x 10~ 5 eV 2 and 
l^mlil = 2.5 x 10~ 3 eV 2 . p is the average matter density along the baseline calculated using the 
Preliminary Earth Reference Model (PREM) [2^ . 



The accuracy of our approximation in calculating the effective mixing angles and effective 
mass-squared differences translates directly into the accuracy in calculating the oscillation 
probabilities. To illustrate this, we calculate the probabilities for the processes — > 
and — > v e for the three cases listed in Table 11111 The energy ranges listed include the 
energies at which IA31J = tt, around where the first oscillation peak occurs. The ranges of 
the matter effect parameter a = log 1 / e (a/|5m| 1 |) for the three cases are roughly —2 ~ — 1, 
— 1.5 ~ —0.5, and —1 ~ 0, so together they cover the range —2 ~ 0. 

Since a < for all three cases, the effective mixing angles are well approximated by 
Eq. ([135)1 . For the effective mass-squared differences, we use Eq. ([146)1 which is applicable 
to a < — 1 for the L = 295 km case, and Eq. () 147)1 which is applicable to a > —1 for the 
L = 2540 km case. 

The L = 1000 km case is a bit problematic since neither Eq. (|146|) nor (|147|) can be used 
throughout the entire range a — — 1.5 ~ —0.5, Eq. ()146j) being applicable only to the low 
energy end, and Eq. ()147jl being applicable only to the high energy end. Using Eqs. ()142j) and 
([143)1 . with Eqs. ([49)1 and ()72j) . which are applicable to all energies, would solve our problem 
and lead to approximate oscillation probabilities that are virtually indistinguishable from 
their exact values. However, we would like to illustrate the power (and limitations) of our 
much simpler expressions in Eqs. ()146j) and ([147)1 . Here, we will use Eq. (1146)1 since the first 
oscillation peak occurs towards the lower end of the energy range. 

In Figs. l23l through 128) we plot the approximate versus the exact oscillation probabilities 
for both the Sm^ > (normal hierarchy) and Sm^i < (inverted hierarchy) cases. The 
vacuum parameters were set to the values listed in Eq. (jl40|) except for the CP violating 
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FIG. 23: Comparison of exact (solid gray line) and approximate (black dashed line) values of P{y^ — ► ^) for the L = 295 km case. The 
approximate value was calculated using Eq. ()135|) for the mixing angles, and Eq. (|146j) for the mass-squared differences. The CP violating 
phase 5 was set to zero. The difference AP = P ap p r ox — -fexact is plotted on the right. 
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FIG. 24: Comparison of exact and approximate values of Piy^ — > v e ) for the L = 295 km case for several different values of the CP violating 
phase 5. The approximate values were calculated using Eq. (|135[) for the mixing angles, and Eq. (|146j) for the mass-squared differences. 
The exact values are given by the solid gray lines, while the approximate values are the black dashed (5 = 0), dotted (5 = tt/2), dot-dashed 
(8 = 7r), and double-dot-dashed (5 = 3n/2) lines. 
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FIG. 25: Comparison of exact (solid gray line) and approximate (black dashed line) values of P{v^ ^u) for the L = 1000 km case. The 
approximate value was calculated using Eq. ()135|) for the mixing angles, and Eq. (|146|) for the mass-squared differences. The CP violating 
phase 5 was set to zero. The difference AP = P ap p r ox — -fexact is plotted on the right. 
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FIG. 26: Comparison of exact and approximate values of -P(f^ — > f e ) for the L = 1000 km case for several different values of the CP violating 
phase 5. The approximate values were calculated using Eq. (|135[) for the mixing angles, and Eq. (|146|) for the mass-squared differences. 
The exact values are given by the solid gray lines, while the approximate values are the black dashed (6 = 0), dotted (5 = 7r/2), dot-dashed 
(5 = vr), and double-dot-dashed (5 = 3ir/2) lines. 
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FIG. 27: Comparison of exact (solid gray line) and approximate (black dashed line) values of P{v^ ^u) for the L = 2540 km case. The 
approximate value was calculated using Eq. ()135|) for the mixing angles, and Eq. (|147j) for the mass-squared differences. The CP violating 
phase 5 was set to zero. 
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FIG. 28: Comparison of exact and approximate values of Piy^ — > f e ) for the L = 2540 km case for several different values of the CP violating 
phase 5. The approximate values were calculated using Eq. (|135[) for the mixing angles, and Eq. (j!47j) for the mass-squared differences. 
The exact values are given by the solid gray lines, while the approximate values are the black dashed (5 = 0), dotted (5 = 7r/2), dot-dashed 
(5 = vr), and double-dot-dashed (5 = 3ir/2) lines. 



phase. For — > i/ e , the four cases 5 = 0, |, 7r, and ^ were plotted. For ^ — > u^, which 
depends only very weakly on S, only the 5 = case is shown. 

The L = 295 km case is shown in Figs. H21 an d E^J As is clear from the figures, our 
approximate values are completely indistinguishable from the exact values, the difference 
being less than a fraction of a percent throughout the energy range considered. For the 
L = 1000 km case shown in Figs. |2H]and|2ni deviations can be seen toward the high energy 
end as expected, but the difference is still well under control. For the L = 2540 km case 
shown in Figs. |2H and |2E1 the difference between the exact and approximate values is again 
less that a percent in the range E = 2 ~ 10 GeV where the approximation is applicable, but 
becomes large below E = 2 GeV. 



V. QUALITATIVE ANALYSIS 

Let us now use our approximation to see whether we can understand various qualitative 
features of the oscillation probabilities we plotted. 

Since the three cases we considered in the previous section cover the range a = —2 ~ 0, 
let us take the midpoint a ~ — 1 and further simplify our expressions to those applicable 
there. The approximations for the effective mixing angles for the neutrinos can be obtained 
from Eqs. (ffti j) . (f7H jl . and (THST) : 







12 



'13 



71 


5mh 


2 ~ 


— 

2a 

6*13 




a 

5m 2 31 


6*23 ) 




5 . 





t?23 

5^5. (175) 

In Fig. I2H1 we plot these approximations against the exact values. Though the approximation 
for 0i3 breaks down for the normal hierarchy case as a — >• 0, these expressions nevertheless 
capture the essential behavior of the effective mixing angles throughout the range — 2 < a < 
0. They and are also numerically accurate for most of this range. In fact, for the inverted 
hierarchy case, the approximations of #13, 023 and 6 are valid for all a. From Eq. (|77jl . the 
approximations for the effective mass-squared differences are 



Ai ~ 5m 



21 °12 1 



Normal Hierarchy Inverted Hierarchy 




FIG. 29: The exact values of #12, #13, O23, and 5 (solid gray lines) plotted as functions of a = 
log 1 / e (a/\5m^i\) against their approximate values (black dashed lines) obtained using Eq. (|175|) . 
The Sm^ > (normal hierarchy) case is shown on the left, and the Sm^ < (inverted hierarchy) 
case is shown on the right. The input parameters are those of Eq. ()14U|) . 
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FIG. 30: The exact values of log 1 / E (<5A2i/|<5m| 1 |) and logi/ e (5Xsi/\Sm1 1 \) (solid gray lines) plotted 
as functions of a = log 1 ^ e (a/\5m'^ 1 1) against their approximate values (black dashed and dot-dashed 
lines) obtained using Eq. ()176|) . The input parameters are those of Eq. (|140j) . 



A 2 ps a + s 2 2 , 

A 3 ~ (Jmfi . (176) 

These expressions are compared against the exact values in Fig. EUJ and we can again 
conclude that they are fairly accurate for most of the range —2 < a < 0, except very near 
the endpoints. 

For the anti-neutrinos, the effective mixing angles when a ~ — 1 are approximated by 

12 ~ ^sin(20 12 ), 
2a 



#13 

013 

1 



Sm 2 31 



6*23 ~ #23 > 

5^5. (177) 

The accuracy of these expressions is shown in Fig. HHJ In contrast to the neutrino case, 
the approximations for # 13 , 9 23 and 5 are applicable to all a for the normal hierarchy case, 
and the approximation for # 13 breaks down as a — > for the inverted hierarchy case. The 
approximations for the effective mass-squared differences are given by 

Ai ~ —a + Sm^ s\ 2 , 
A 2 ps 5m 2 21 c\ 2 , 

A 3 ~ 5m 2 31 , (178) 

with the accuracy shown in Fig. (JH2*j) . As in the neutrino case, we can conclude that these ap- 
proximations capture the essential behavior of the effective mixing angles and mass squared 
differences throughout the range —2 < a < 0, and are also numerically accurate except near 
the endpoints. 

Let us now apply these approximations to the oscillation probabilities Piy^ — > v^) and 
P{yn — ¥ v e) m matter, and their anti-neutrino counterparts. First, recall from Eqs. I}2()j) and 
(12 lj) that these probabilities in vacuum are given by 

P{v, - v.) = 1 - sin 2 (2^ atm ) sin 2 ^ 31 ~ + 0(A 2 X ) , 

P{v^v e ) = 4sin 2 £ 13 si^^jl-^sin^A^jsin 2 ^ ~ ^A 21 ^ + 

(179) 
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FIG. 31: The exact values of #12, #13, O23, and 5 (solid gray lines) plotted as functions of a = 
log 1 / e (a/\5m^i\) against their approximate values (black dashed lines) obtained using Eq. (|177|) . 
The Sm^ > (normal hierarchy) case is shown on the left, and the Sm^ < (inverted hierarchy) 
case is shown on the right. The input parameters are those of Eq. ()14U|) . 
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FIG. 32: The exact values of log 1 / E (<5A2i/|<5m| 1 1) and logi/ed^Asij/IJm^l) (solid gray lines) plotted 
as functions of a = log 1 ^ e (a/\5m'^ 1 \) against their approximate values (black dashed and dot-dashed 
lines), obtained using Eq. (|178|) . The input parameters are those of Eq. (|140|) . 



where 



sin atm = S23C13 , 

A = I sin(2# 12 ) sin(2# 13 ) sin(20 atm ) y/l - tan 2 13 tan 2 atm 

o 

= - x 0(1) x 20(e) x 0(1) x y/l - 0(e 2 ) 

= > 

sm 013 sm atm 2 
K w = c i2 - 608(2012) tan 2 0i3 tan 2 atm - I — — — J cos 5 

\C0S 2 013 COS atm / 

= - 0(1) -20(e 3 ) -0(e) cos5 
= - 0(1) -0(e) cos 5 , 

K A*e = s 12 ~ B COS S 

= iocij-iocOcos*. (180) 

The oscillation probabilities for the anti-neutrinos are obtained by flipping the sign of sin 5. 

Note that the first oscillation peak occurs at a distance /energy of IA31I ~ it. Since 
A 2 i/|A 3 i| = Sm^/lSm^l = e 2 , the O(A^) terms in Eq. (ITTHll are of 0(vr 2 e 4 ) which justifies 
our dropping them at those distance/energies. Note also that the coefficient of cos 5 in K^ e 
is B = ^(e" 1 ), which when multiplied by A21 is of 0(e). The exact same product of 
parameters, -BA 2 i = 0(e), appears in the coefficient of sin 5 in the oscillation envelope of 
P{ v n ~^ v e)- Therefore, a measurement of P{v^ — > u e ) can, in principle, constrain cos S 
from the position of the peak, and sin 5 from the height of the peak, provided it is accurate 
enough to discern these 0(e) corrections. In contrast, the coefficient of cos 5 in is 0(e), 
rendering the cos 5 term in k m/j A 2 i negligible, and we conclude that Piv^ — * Vp) is insensitive 
to 5. 

Of course, actual long-baseline experiments can only measure oscillation probabilities in 
matter. The effective A 3 i's in matter for the neutrinos and anti- neutrinos are 



A 3 i « (A 3 i-A 2 ic 2 2 ) 

respectively, so they are both the same order as A 3i . The effective A 2 i's, on the other hand, 



A 3 i « (A31 + - A 2 is 2 2 ) , (181) 



arc 



A 2 i « (^L-A 2 i cos 20i 2 ) « = (V2G F N e )L 



A 2 i « + A 21 cos2# 12 ) « = (y/2G F N e )L , (182) 

so they are enhanced by a factor of a/Sm^ relative to A21. Therefore, at the first oscillation 
peak where |A 31 | m it, we can expect A 2 i and A 2 i to be of order ita/\5m\ x \. 

Note that the value of A 21 ps A 21 « (\^2GpN e )L does not depend on the energy E 1 . It is 
determined solely by the baseline length L, once the matter density p is fixed. (N e = 
where Na is the Avogadro number.) For the three examples we considered in the previous 
section, we find 



{y/2G F N e )L = < 



0.15 (p = 2.6 g/cm 3 ,L = 295 km) , 
0.5 = 2.7 g/cm 3 ,L = 1000km) , (183) 
1.7 (p = 3.4 g/cm 3 ,L = 2540km) . 

Since we would like to use the analog of Eq. (| 179(1 to analyze the oscillation probabilities in 
matter, we would like maintain the condition 

A 21 , A 21 < 1 , (184) 

so that an expansion in A 2 i or A 2 i is justified. The L = 2540 km case is clearly problematic 
and must be treated separately. We will therefore first restrict our attention to the cases 
in which the |A 31 | = it condition occurs in the region —2 < a < —1, so that A 21 A 21 < 
0(ne) = 0.47 ~ 0.75. 

Even with this restriction, A 21 and A 21 can still be enhanced considerably when |A 31 | = 7r 
occurs at a ~ —1. One may naively anticipate that this enhancement will enhance the 
coefficients of sin S and cos 5 in Piv^ — > p e ), thereby facilitate the detection of 5. At the 
same time, it could also enhance the 0(A 21 ) and 0(A 21 ) terms in the oscillation probabilities 
and invalidate their complete neglect. However, it turns out that these are not the case. 

Let us first look at the and survival probabilities in matter which are obtained 
by replacing all the quantities in the vacuum probability with their tilded and anti-tilded 
counterparts: 

= 1 - sin (20 atm ) sin I 2 

l^iH^I 2 (l+ ^ Ulx f 2 eosA 3 i 1 A 21 



1 - U 



fj3 1 




smA 3 A A 21 + 0(A 21 ) , (185) 



fj3 1 



sinA 3 i } A 21 + 0(A 21 ) . (186) 



We have kept terms up to A 21 and A 21 explicitly to evaluate their sizes. Using the approx- 
imations of Eqs. ()175j) and (j!77)l . the effective MNS matrix elements that appear in these 
expressions can be evaluated to be 

~ ~ i 

\U^\ 2 = \s 12 c 23 + ci 2 si3S2 3 e J,5 | 2 = -0(1) , 

1 2 |~ ~ ~ ~ ~ „iS\2 



\Ufa\ = |Cl 2 C 2 3 - 5125135236* | « 0(e 

2 



\U, 3 \ 2 = (ci 3 5 23 ) 2 = lo(l) . (187) 



for the neutrinos, and 



\Upi\ = \s 12 C23 + ci2Si 3 S2 3 e I « 0(e), 
~ ^ 1 

I T r |2 , ^ ^ ^ ^ ^ 7 A i 2 

|C/ M 2| = |Cl 2 C 2 3 - Si 2 Si3S2 3 e | = T^V 1 ) > 

IM 2 = (c 13 3 23 ) 2 = \o{\) , (188) 

for the anti-neutrinos. This shows that the A 21 , A 21 and A 21 , A 21 terms are suppressed by 
I U n2 1 2 = 0(e 2 ) for the neutrinos, and by |^i| 2 = 0(e 2 ) for the anti-neutrinos, cancelling out 
the enhancements of A 2i , A 21 over A 21 , and reducing these higher order terms to 0(7r 2 £ 4 ), 
which is the same order as the 0(A 21 ) terms that were neglected for the vacuum case. The 
fact that cos A 31 w — 1 and sin A 31 pa near the first oscillation peak also helps in suppressing 
these terms. Therefore, to a very good approximation, we can use the expressions 

P(vn -»• v^) = 1 - sin 2 (26» atm ) sin 2 

- ^) = 1 - sin 2 (2^ atm ) sin 2 I " 31 ' W ~ Z1 ) . 




Looking at the remaining parameters in these expressions, we find 
sin^atm = C13S23 = C13S23 [l + 0(e 3 )] « sm6» atm , 



A = ^sin(2^i2)sin(2^i3)sin(2^ atm )\/l - tan 2 # 13 tan 2 ato 

o 

1 ( Smli 



A \ a J \} Sm^J ° 13 

2 ~ ~ ~ / 2y~l 
*W = ^12 _ cos(26 l 1 2) tan 2 13 tan 2 atin - — J cos 5 

\COS 2 0i 3 COS 2 8 E 

= 0(s 2 ) -0(e 2 ) cos 8 , (190) 
for the neutrinos, and 

sin0 atm = C13S23 = C13S23 [l + 0(e 3 )] « sin atm , 



A « -sin(2^ 1 2)sin(2e 13 )sin(2e atm )\/l - tan 2 13 tan 2 



' atin 



1 / 5m 2 ,, \ { a \ . 



k w = cos 2 i2 - cos(26 l i2) tan 2 13 tan 2 atm - - — J cos 5 

\ cos 2 0i3 cos 2 atm / 
= 0(1) -0(e 2 )cos5 , (191) 

for the anti-neutrinos. (Recall that sin(26 l i2) = 1 — 20(e 2 ) and sin(26 l atm ) = 1 — \0(e 2 ).) 
This shows that the coefficient of cos 5 in (iS^) is suppressed by one power of e relative to 
K MAt . Consequently, the and z/ M survival probabilities remain insensitive to the CP violating 
phase 5, despite the enhancements of A21, A21 over A21. Looking at the arguments of the 
sines more carefully, we find 

A 3 i - « w A 2 i = (A 3 i - c 2 2 A 2 i) + 0(e 3 A 31 ) , 

A 3 i - K wt A 2i « (A 3 i - c 2 2 A 2 i) - c 2 12 (wjjjL ~ A 2i cos20 12 ) 

= (A 3 i-c 2 2 A2i)+0(£ 3 A 3 i), 

A 3 i - /? MM A 2 i « (A 31 + - A 2 is 2 2 ) - cj 2 (^L + A 2 i cos20 12 ) 

= (a 3 i + - A 2 is 2 2 ) - (^L + A 2 i cos 20 12 ) + 0(s 3 A 31 ) 

= (A 3 i-c 2 2 A2i)+0(£ 3 A 3 i). (192) 

Since sin(26 l a t m ) ~ sin(20 atm ) ~ sin(20 atni ), we conclude 

P(v» -> ^) « P(^ ^) = P(^ P„) ~ P(^ M -> ^) • (193) 



That is, the and survival probabilities are insensitive to matter effects and their 
values in matter are the same as their values in vacuum. In Fig. EH1 we compare the exact 
numerical values of these probabilities calculated for our example parameter set Eq. ()140J) 
with Sm^ > 0. As can be seen, the differences among these probabilities are extremely 
small, and our approximation has allowed us to understand this analytically. 

Next, let us consider the — > v e and — > v e oscillation probabilities. They are 



4 sin 2 0i3 sin 2 E 



(5 sin5)A 2 i|sin 2 



+ - < (B sin 5) 2 - 2« Me (l - n^ e ) sin 2 — 



^31 — ^/ie^21 



(5 sin 5) (1 - sin A 31 > A 21 



~ i 3(5 sin 5)%l e + 2(5 sin 5) (1 - 3«Jj sin 2 ^ - K, e (l - sin A 3 i j> A 



3 
21 



+ 0(A 



21- 



(194) 



4 sin 2 0i3 sin 2 



' atm 



| A 3 i — K Me A 2 l 

1 + (5 sm o)A 2 i f sin 1 



+- <( (5sin<5) 2 - 2/^(1 - k^) sin 2 ^- + (5 sin 5) (1 - 2/? Me ) sin A 31 A 21 



3(5 sin £)^ e + 2(5 sin 6)(1 - 3«£j sin 2 + £ Me (l - «£j sin A 3i 



+ 0(A 



21. 



(195) 



From Eq. (|19U|) we obtain 



5 



'7<C:' 



.4 



1 f dm^ 



sin 2 0i3 sin 2 8 atm 



1 



2 0(1) 



~2 1 

s 12 — 5 cos 5 ~ 1 — - 0(1) cos 5 



and from Eq. (|191|) we obtain 
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FIG. 33: Comparison of the and survival probabilities in vacuum and in matter. On the left, Piy^ — ► v^) = P{pa £v) is the solid 
gray line, P(vu — ► v u) is the dashed black line, and Pipp — > Va) is the dotted black line. On the right, the difference AP(z/„ — ► Un) = 
P{V[i — > Vp) — Piy^i — ► Vn) is the dashed black line, and the difference AP(^ — > U^) = P{pp — > Vp) — Pipp, —> v^) is the dotted blackline. 
The input parameters are those listed in Eq. ()14U|) with Sm^ > 0. 



K 



f 12 -Bcos8 = ]o(e 2 ) - ^O(l)coso" « --0(l)cos8 . 



i-sj 2 2 3 ^-^ 3 

Therefore, the expressions inside the curly brackets in the A 21 , A 21 and A 21 , A 21 terms are 
roughly of order one, with factors of 1/4 and 1/12 in front suppressing them. And since the 
entire expression is multiplied by sin 2 #13 = 0(e 2 ), sin 2 13 = 0(e 2 ) from outside the square 
brackets, these terms are of 0(7r 2 e 4 ) and are negligible. 

Therefore, to a good approximation, we can use the expressions 



4 sin 2 0i3 sin 2 atm \ 1 - (B sin 8) A 2 i f sin 



: 4 sin 2 0i3 sin 2 atm { 1 + (B sin 5) A 21 } 
where the arguments of the sine functions are given by 



sin 




(196) 



A 



31 



fvj A 

KupLA 



A31- 

(A 31 

A 3 i 

A 3 i 



Si 2 — 5 COS 8 ) A21 



1 / 5m 21 



! V a0i 3 
20i, Vom 2 



8m 2 31 



~ ^21 cos 20 i2 



sign(5m 2 1 ) 



cos 5 + 0(e' 

r-2 



\8m 2 



311 



20 



13 



cos 8 J + 0{e" 



A31 — K Me A 2 i 



A 



31 



12 

a 

2£" 



I? cos 8 ) A 2 i 



A31 + — L - s 2 2 A 2 i 



a0i 3 



1 + 



A31 



8m 2 31 



20 



13 



pf) C os8 + 0(e 2 ' 



4) 



\8m 



2 1 

31 1 



sign(o~m 

while the amplitudes are given by 



20 



cos 8 



13 



0(e* 



cos 8 



— L + A 2 i cos 20i2 



(197) 



20i3 



20 2 3 



2a 

2a 
omli 



- (fisin5)A 2 i| 


1 - < 


fl /5m 2 , 




1 2 Va0i3 




^m 2 A A 


20i3 ( 


v om 2 J 



1 - 



8m 2 31 



sin 8 



A 3 i sin 5 + 0(e" 



/ a 
\2~E 



L — A 2 i cos 20i2 



2 #13 



2a 
8m 2 31 



20 



A 3 i| sin 8 + 0(e* 



13 



4sin^i 3 sin^ atm <! 

2a 



26> 2 3 



26> 2 3 



2 ^13 



1 - 
1 - 
1 - 



5m 2 31 
2a 

2a 



(Bsm5)A 2 i 
1 
2 



aO 



13 



1 / 5m 2 



A 2 l cos 26> 12 



20 



21 



A31 sin 5 + 0(e 2 



20 



A31I sino" + Ofe 5 



13 



'1981 



From these simple expressions, we can discern a few facts about v e (i> e ) appearance experi- 
ments. 

Roughly speaking, the positions of the oscillation peaks will provide information on the 
combination 



sign(om 



=F 



a 



\5ml 



20, 



cos 8 



(199) 



'311 

while the heights of the oscillation peaks will provide information on the combination 



20\ 3 



1 ± < sign (8m 



2a 



A31I sin 8 



(200) 



\8m 2 31 \ 20 13 ' 

where the upper signs are for the neutrinos, and the lower signs are for the anti- neutrinos. 
If the oscillation peaks occur in an energy region in which 

-2 



\8m 2 



< 



£ 



20i 



(201) 



'311 ^13 

(a ~ —2) then the a/loYnf-J terms in these expressions can be neglected. Then, measuring 
the height of the peak (or just the total v e or v e flux using a narrow band beam) will allow us 



27 



to constrain sin 8, provided that #13 is well known from future reactor experiments 
If 9 13 is not well-known, then measuring the peak heights for both the neutrino and anti- 
neutrino will allow us to constrain both # 13 and sin 8. If one also measures the position of 
the peaks, either using a wide band beam or by changing the beam energy, then one can 
also extract information on the product 



sign(o"m3 1 ) cos 8 



(202) 



but neither sign(o"mg 1 ) nor the sign of cos 8 can be uniquely determined, even if both neutrino 
and anti-neutrino beams are used 



These features are clearly visible in Fig. |2U which shows the probabilities to be probed by 
the T2K experiment Q|. In phase 2 of T2K, the oscillation event rates for both neutrinos 
and anti-neutrinos are to be measured, the difference from which we can extract sin 8. 
However, this does not provide any information on cos 5. It was proposed in Ref. H to 
measure the event rates using several beams of different energy, and thereby obtain some 
information on the peak position of the oscillation spectrum, but the sign of cos 5 cannot be 



uniquely determined unless the sign of Sm^ is known 31 1. (These difficulties can be best 
seen visually by utilizing the Minakata-Nunokawa plot 32j.) 

On the other hand, if the oscillation peaks occur in an energy region in which 

a e 2 

¥<\ w Ws ' (203) 

(a ~ —1) then the measurement of the peak height by itself may not be able to determine 
either sign^m^) or sin 5, even if #13 were accurately known. In particular, there will be a 
degeneracy between the two cases in which sign(5mg 1 ) and sin S are both positive, and both 
negative, as can be clearly seen in Fig. I2E1 Measuring peak heights for both the neutrino and 
anti-neutrino will not help use here since they both depend on the same linear combination 
of a/drn?^ and sin 5, though it will help us in determining # 13 . If one also measures the 
position of the peak, then the degeneracy in sign^m^) and cos 5 which existed for the 
previous case can be lifted, except when <5 ^ for the neutrino case, and 5 ~ n for the 
anti-neutrino case. Due to these shortcomings in performing a single experiment at either 
a ~ —1 or a ~ —2, various scenarios have been suggested which utilize two detectors set 
up .t different baseUne ien^s QQQQ. 

For the L = 2540 km case, an expansion in A 2 i is no longer permissible. The approxi- 
mation we used above for 13 also breaks down as a — > 0. We can nevertheless simplify the 
expressions for the oscillation probabilities and understand their behavior analytically. As 
an example, consider the oscillation probabiliity Piv^ — > v e )-, the full expression of which is 

= 4|^ 2 | 2 |^ 2 | 2 sin 2 ^i + 4|^ 3 | 2 |^3| 2 sin 2A31 



2 

L 

2 " 

+4 J M I sin 2 sin A 31 - sin 2 =£i sin A 21 j . (204) 



+2 ^{U^U^U^Ua) \ 4 sin 2 ^ sin 2 ^ + sin A 21 sin A 31 
.2 ^21 . 7 . 2 ^31 



Using the approximations in Eq. (J175j) . except the one for #i 3 , we find 



\Ufi2\ 2 \U e 2\ 2 ~ S 23 S 13 C 13 — 2 A COS 5 , 

I tt |2|rr 1 2 2 ^2 

|^Ai3| I^e3| ~ S 23 S 13 C 13 > 

e2 ) ~ -S^S 13 C 13 + ACOS5 , 

J^ej « 2 sin 5 , (205) 



where 



5m 2 



A « ( ) sin(2^ 12 )s 2 3C 2 3S 13 cJ3 . (206) 



Terms of order (o~m 21 /a) 2 and higher have been neglected. Substituting into Eq. (|204|h we 
obtain 
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■ foa \ ■ A 32 ~ 
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'A 31 




sm cos 




V a , 


' 2 \ 


, 2 





(207) 



where we have used sin(26'i 2 ) = 1 — 20(e), and sin(20 23 ) = 1 — \0(e 2 ). Now, the a- 
dependence of #13 is different depending on the sign of Sm^. If Sm^ > (normal hierarchy), 
#13 increases monotonically from #13 toward | as a increases with energy, and passes through 
j around a « joVnl-J. If Sm^ < (inverted hierarchy), # 13 decreases monotonically from 
813 toward as a increase with energy, and is about #13 /2 around a m |5m 31 |. (cf. Table ITVj) 
Therefore, if 5m 31 > 0, then the coefficient of the first term in the brackets of Eq. ()207|) is 
maximized around a ~ l^m^-J. If the oscillation peak where A 32 m n matches that energy, 
one can expect a maximum oscillation probability as large as |. (Note that this maximum 
probability is determined by # 23 and is independent of the value of 6*13 in vacuum.) On the 
other hand, if Sm^ < 0, then the same coefficient is suppressed to 0(e), and the oscillation 
probability with be suppressed by a factor of 0(e 2 ). This difference is evident in Fig. [2B1 
(In both cases, the second term in the brackets interferes with the first term giving the 
probability a weak (^-dependence.) Due to this clear difference, measuring the v e appearance 



probability at the first oscillation peak at a baseline length o 



as an unambiguous method to determine the sign of 5m 31 



L = 2540 km has been proposed 
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TABLE IV: The dependence of the effective mixing angles on a/l^m 2 ^. 

Comparing the L = 295 km, L = 1000 km, and L = 2540 km cases considered above, we 
can discern a generic trend that the oscillation probability Piv^ — > v e ) (or Piv^ — > v e )) is 
more sensitive to the CP violating phase 5 at lower energies (shorter baselines), and more 
sensitive to the mass hierarchy (sign of Sm 3l ) at higher energies (longer baselines). This 
can be understood as follows. First, the angle #13, which has a different energy dependence 
depending on the sign of Sm 31 , enters Pivp, — > v e ) dominantly in the combination 

sin 9 13 sin # atm = sin 9 13 cos 9 13 sin 6 23 
= -sin(26>i 3 ) sin ^23 

« ^sin(2# 13 )sin# 23 . (208) 

Therefore, sin(26 , 13 ) determines the sensitivity of Piv^ — > v e ) on the behavior of # 13 . On 
the other hand, the size of CP violation is governed by the Jarskog invariant, 

J = Asm 5, (209) 

where 

~j _ ~ ~ ~ ^2^ ^ 

A — Si 2 Ci 2 Si 3 C 13 S2 3 C2 3 

1 ~ ~ ~2 

= - sin(26»i 2 ) sin(26> 23 )si 3 (l - s 13 ) 

« i sin(2# 12 ) sin(2# 23 )£ 13 (l - s' 3 ) . (210) 
In our convention where the mixing angles are in the first quadrant, A is bounded by 




FIG. 34: The dependence of sin(2#i3), sin(2#i3) (solid gray lines), sin(2#i2), sin(26>i2) (dashed 
black lines), £13(1 — s 13 ), 813(1 — s 13 ) (dotted black lines), and A/A max , A/A maiyi (solid black lines) 
on a = log 1 / e (a/|5?n| 1 |). The input parameters were those listed in Eq. (|14()|) . 



In Fig. EU we plot the a-dependence of sin(26 l i3), sm(28 12 ), £13(1 — s 31 ), and A/A mSLX , and 
their anti-neutrino counterparts, for the input parameters of Eq. (jl4(J|) . It is clear from the 
figure that the differences in sin(2# 13 ) and sin(2# 13 ) between the Sm^ > and Sm^ < 
cases is most prominent at a ~ where there exists a large peak in one case which is absent 
in the other. This is due to the fact #13 (#13) crosses n/A at a m {Sm^] when Sm^ > 
(Sm^i < 0). (cf. Table ITVh Therefore, experiments that are performed closer to a = are 
more sensitive to sign^m^). On the other hand, the a-dependence of A (A) is dominated 
by that of sin(26 l 1 2) (sin(2# 12 ) which starts out as a number of 0(1) at low a, but drops 
off quickly toward zero as a is increased from —2 to —1. Therefore, experiments that are 
performed at a < —2 are more sensitive to S than those performed at a > — 1. 



VI. SUMMARY 



In this paper, we considered the matter effect on neutrino oscillations, and derived simple 
analytical approximations to the effective mixing angles and effective mass-squared differ- 
ences in constant density matter. Our results are summarized in Table IVl These expressions 
can be utilized in calculating, analyzing, and understanding the behavior of the oscillation 
probabilities in LBL neutrino oscillation experiments. 

The formalism developed in this paper can be further extended to incorporate matter 
effects due to the violation of universality in neutral current interactions, or additional 
matter effects due to new interactions. This will be presented in a subsequent paper [3^ . 
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TABLE V: The approximate formulae for the effective mixing angles and effective mass-squared differences derived in this paper. 
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